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Summary 



In this thesis, we consider semi-algebraic sets over a real closed field R 
defined by quadratic polynomials. Semi-algebraic sets of R'^ are defined as 
the smallest family of sets in R'^ that contains the algebraic sets as well as the 
sets defined by polynomial inequalities, and which is also closed under the 
boolean operations (complementation, finite unions and finite intersections). 
We prove the following new bounds on the topological complexity of semi- 
algebraic sets over a real closed field R defined by quadratic polynomials, in 
terms of the parameters of the system of polynomials defining them, which 
improve the known results. 

(1) Let 5 C R^ be defined by Pi > 0, . . . , > with 

Pi £ R[Xi, . . .,Xk], m < k, and deg(Pi) < 2, for 1 < i < m. We 
prove that bi{S) < § • (^)"' + A;, < i < A; - 1. 

(2) Let V = {Pi,..., Pm} C R[Yi, ...,Ye,Xi,.. .,Xk], 

with degy(Pi) < 2, degx(Pi) < d, I < i < m. Let S" C R^+'' be a 
semi-algebaic set, defined by a Boolean formula without negations, 
whose atoms are of the form, P > 0, P < 0, P G 7^. 
Let vr : R^"'"'' — > R'"' be the projection on the last k co-ordinates. 
We prove that the number of stable homotopy types amongst the 
fibers 7r-i(x) n 5 is bounded by 
We conclude the thesis with presenting two new algorithms along with their 
implementations. The first algorithm computes the number of connected 
components and the first Betti number of a semi-algebraic set defined by 
compact objects in M'^ which are simply connected. This algorithm improves 
the well- know method using a triangulation of the semi- algebraic set. More- 
over, the algorithm has been efficiently implemented which was not possible 
before. The second algorithm computes efficiently the real intersection of 
three quadratic surfaces in M'^ using a semi-numerical approach. 
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CHAPTER 1 



Introduction 

1. Real Algebraic Geometry 

In classical algebraic geometry, the main objects of interest are complex 
algebraic sets, i.e. the zero set of a finite family of polynomials over the 
field C of complex numbers, meaning the set of all points that simultane- 
ously satisfy one or more polynomial equations. But in many applications 
in computer-aided geometric design, computational geometry, robotics or 
computer graphics one is interested in the solutions over the field M of real 
numbers. Moreover, they also deal with the real solutions of finite systems of 
inequalities which are the main objects of real algebraic geometry. Unfortu- 
nately, real algebraic sets have a very different behavior than their complex 
counterparts. For example, an irreducible algebraic subset of C'^ having 
complex dimension n, considered as an algebraic subset of M."^^ is connected, 
not bounded (unless it is a point) and has local real dimension 2n at every 
point (see, for instance, |27j ). But this is no longer true for real algebraic 
sets (see Example 12. 38p . 

In 1926, Emil Artin and Otto Schreier [Tl [6] introduced the notion of 
a real closed field. Artin [5|, [6] used this new theory for solving the 17th 
problem of Hilbert which asks whether a polynomial which is nonnegative 
on M" is a sum of squares of rational functions. A real closed field R is an 
ordered field whose positive cone is the set of squares R^^^ and such that 
every polynomial in R[X] of odd degree has a root in R. Notice that real 
closed fields need not be complete nor archimedean (see Chapter II. 2p . 

In this thesis, we consider semi-algebraic sets over a real closed field R 
defined by quadratic polynomials in k variables. Semi-algebraic sets of R'^ 
are defined as the smallest family of sets in R'^ that contains the algebraic 
sets as well as the sets defined by polynomial inequalities, and which is 
also closed under the boolean operations (complementation, finite unions 
and finite intersections). Furthermore, unlike algebraic sets (over R), the 
projection of a semi-algebraic set is again semi-algebraic, this was proved by 
Tarski [75] and Seidenberg [70] . 

It is worthwhile to mention that in many applications in computer-aided 
geometric design or computational geometry one deals with arrangements 
of many geometric objects having a similar simple description j48j . For 
instance, each object is a semi-algebraic set defined by few polynomials of 
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2 1. INTRODUCTION 

fixed degree. Thus, understanding the properties of semi-algebraic sets and 
designing algorithms are important topics in real algebraic geometry. 

The class of semi-algebraic set defined by quadratic polynomials is of 
particular interest for several reasons. First, any semi-algebraic set can be 
defined by (quantified) formulas involving only quadratic polynomials (at 
the cost of increasing the number of variables and the size of the formula). 
Secondly, they are distinguished from arbitrary semi-algebraic sets since one 
can obtain better results from an algorithmic standpoint, as well as from 
the point of view of topological complexity (as we will see later). Moreover, 
they can be much more complicated topologically than semi-algebraic sets 
defined by only linear polynomials. Thirdly, quadratic surfaces are widely 
used in computer-aided geometric design, computational geometry j69j and 
computer graphics as well as in robotics ( |68j ) and computational physics 

(mm)- 

One basic ingredient in most algorithms for computing topological prop- 
erties of semi-algebraic sets is an algorithm due to Collins [33], called cylin- 
drical decomposition (see Chapter II. 4p which decomposes a given semi- 
algebraic set into topological balls. Cylindrical decomposition can be used 
to compute a semi-algebraic triangulation of a semi-algebraic set (see Chap- 
ter [T3]), and from this triangulation one can compute the homology groups, 
Betti numbers, et cetera. One disadvantage of the cylindrical decomposition 
is that it uses iterated projections (reducing the dimension by one in each 
step) and the number of polynomials (as well as the degrees) is squared in 
each step of the process. Thus, the complexity of performing cylindrical de- 
composition is double exponential in the number of variables which makes it 
impractical in most cases for computing topological information. Neverthe- 
less, we will see in Chapters 11.4.21 and [51 that it can be used quite efficiently 
for several important problems in low dimensions. 

2. Betti numbers 

Important topological invariants of a semi-algebraic sets are the Betti 
numbers bi (see Chapter 12. II for a precise definition) which, roughly speaking, 
measure the number of i-dimensional holes of a semi-algebriac set. The zero- 
th Betti number bo is the number of connected components. 

The initial result on bounding the Betti numbers of semi-algebraic sets 
defined by polynomial inequalities was proved independently by Oleinik and 
Petrovskii [65], Thorn |76| and Milnor [63] . They proved (see Theorem l2.34p 
that the sum of the Betti numbers of a semi-algebraic set in R'^ defined by 
m polynomial inequalities of degree at most d has a bound of the form 
0{md)^ . Notice that this bound is exponential in k and this exponential 
dependence is unavoidable (see Example 12.35^ . Recently, the above bound 
was extended to more general classes of semi-algebraic sets. For example, 
Basu [llj improved the bound of the individual Betti numbers of P-closed 
semi-algebraic sets (which are defined by a Boolean formula with atoms 
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of the form P = 0, P < or P > 0, where P (z V), while Gabrielov and 
Vorobjov |44j extended the above bound to any P-semi-algebraic set (which 
is defined by a Boolean formula with atoms of the form P = 0, P < or 
P > 0, where P G V). They proved a bound of 0{m?d)^. Moreover, 
Basu, Pollack and Roy |19j proved a similar bound for the individual Betti 
numbers of the realizations of sign conditions. 

However, it turns out that for a semi-algebraic set 5 C R'' defined by 
m quadratic inequalities, it is possible to obtain upper bounds on the sum 
of Betti numbers of S which are polynomial in k and exponential only in 
m. The first such result was proved by Barvinok [9] who proved a bound of 
)^0{m) ^ggg Theorem l2.36p . The exponential dependence on m is unavoidable 
as already remarked by Barvinok, but the implied constant (which is at least 
two) in the exponent of Barvinok's bound is not optimal. 

Using Barvinok's result, as well as inequalities derived from the Mayer- 
Vietoris sequence (see Chapter 12. 2p . Basu proved a polynomial bound 
(polynomial both in k and m) on the top few Betti numbers of a set defined 
by quadratic inequalities (see Theorem l2.37p . Very recently, Basu, Pasechnik 
and Roy |18j extended these bounds to arbitrary 7^-closed (not just basic 
closed) semi-algebraic sets defined in terms of quadratic inequalities. 

Apart from their intrinsic mathematical interest, for example in dis- 
tinguishing the semi-algebraic sets defined by quadratic inequalities from 
general semi-algebraic sets, the bounds proved by Barvinok and Basu re- 
spectively have motivated recent work on designing polynomial time algo- 
rithms for computing topological invariants of semi-algebraic sets defined by 
quadratic inequalities. For instance, Grigoriev and Pasechnik |47j presented 
a polynomial time algorithm (in k) for computing sampling points meeting 
each connected component of a real algebraic set defined over a quadratic 
map. Their result improves a result of Barvinok [8] about the the feasibil- 
ity of systems of real quadratic equations. Basu |14|, I13j gave polynomial 
time algorithms for computing the Euler characteristic and the higher Betti 
numbers of semi-algebraic sets defined by quadratic inequalities. Further- 
more, Basu and Zell [23j gave a polynomial time algorithm for computing 
the lower Betti numbers of projections defined by such semi-algebraic sets. 
For details, we refer the reader to the papers mentioned above. 

Traditionally an important goal in algorithmic semi-algebraic geometry 
has been to design algorithms for computing topological invariants of semi- 
algebraic sets, whose worst-case complexity matches the best upper bounds 
known for the quantity being computed. It is thus of interest to tighten the 
bounds on the Betti numbers of semi-algebraic sets defined by quadratic in- 
equalities, as it has been done recently in the case of general semi-algebraic 
sets (see for example [44^ 1111, I19^ 118] ). Notice that the problem of com- 
puting the Betti numbers of semi-algebraic sets in single exponential time 
is considered to be a very important open problem in algorithmic semi- 
algebraic geometry. Recent progress has been made in several special cases 
(see [2Tl[T2t[T4]). 
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In another direction, the bounds of the Betti numbers are used to pro- 
duce lower bounds for complexity decision problems. For instance, Steele 
and Yao |74j recognized that the bounds for the sum of the Betti numbers 
can be applied to obtain non-trivial lower bounds in terms of the number of 
connected components for the model of algebraic decision trees. This was 
extended to algebraic computation trees by Ben-Or |25j . 



3. Homotopy Types 

A fundamental theorem in semi-algebraic geometry is Hardt's Theorem 
(see Theorem l2.15p which is a corollary of the existence of the cylindrical de- 
composition. For a projection map tt : R^^'^ — > R'^ on the last k co-ordinates 
and semi-algebraic subset S of R'^, it implies that there is a semi-algebraic 
partition of R'^, {Tj}jg/, such that for each i G / and any point y S Tj, the 
pre-image iT~^{Ti)r\S is semi-algebraically homeomorphic to (7r~^(y)nS') xTj 
by a fiber preserving homeomorphism. In particular, for each i £ I, all fibers 
7r~^(y) n S, y € Tj, are semi-algebraically homeomorphic. Unfortunately, 
the cylindrical decomposition algorithm implies a double exponential (in k 
and i) upper bound on the cardinality of / and, hence, on the number of 
homeomorphism types of the fibers of the map vrl^. No better bounds than 
the double exponential bound are known, even though it seems reasonable 
to conjecture a single exponential upper bound on the number of homeo- 
morphism types of the fibers of the map vr^. 

Basu and Vorobjov |22j considered the weaker problem of bounding the 
number of distinct homotopy types, occurring amongst the set of all fibers 
of -7r|5, and a single exponential upper bound was proved on the number of 
homotopy types of such fibers (see Theorem 12. 42p . They proved in the same 
paper a similar result for semi-Pfaffian sets as well, and Basu |14] extended 
it to arbitrary o-minimal structures. Both these bounds on the number of 
homotopy types are exponential in i as well as k. As already pointed out in 
|22j . in this generality the single exponential dependence on i is unavoidable 
(see Example [233]). 

Since sets defined by quadratic equalities and inequalities are the sim- 
plest class of topologically non-trivial semi- algebraic sets, the problem of 
classifying such sets topologically has attracted the attention of many re- 
searchers. Motivated by problems related to stability of maps. Wall [79] 
considered the special case of real algebraic sets defined by two simultane- 
ously diagonalizable quadratic forms in i variables. He obtained a full topo- 
logical classification of such varieties making use of Gale diagrams (from the 
theory of convex polytopes). To be more precise, letting 
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i=l 

and 

5 = {(y,x) gR3^| II y 11=1, Qi(y,x) =Q2(y,x) =0}, 

Wall obtains as a consequence of his classification theorem, that the number 
of different topological types of fibers 7r~^(x)n5' is bounded by 2^~^. Similar 
results were also obtained by Lopez [59j using different techniques. Much 
more recently Briand |28j has obtained explicit characterization of the iso- 
topy classes of real varieties defined by two general conies in two dimensional 
real projective space Pp^ in terms of the coefficients of the polynomials. His 
method also gives a decision algorithm for testing whether two such given 
varieties are isotopic. 

In another direction Agrachev [l] studied the topology of semi-algebraic 
sets defined by quadratic inequalities, and he defined a certain spectral se- 
quence converging to the homology groups of such sets. We will give a 
parametrized version of Agrachev's construction in Chapter [3] which is due 
to Basu. 

In view of the topological simplicity of semi-algebraic sets defined by few 
quadratic inequalities as opposed to general semi-algebraic sets, one might 
expect a much tighter bound on the number of topological types compared 
to the general case. However one should be cautious, since a tight bound on 
the Betti numbers of a class of semi-algebraic sets does not automatically 
imply a similar bound on the number of topological or even homotopy types 
occurring in that class. We refer the reader to |24j for an explicit example of 
the large number of possible homotopy types amongst finite cell complexes 
having very small Betti numbers. 

4. Arrangements 

Arrangements of geometric objects in fixed dimensional Euclidean space 
are fundamental objects in computational geometry and computer-aided 
geometric design (for instance, see [48]). As already mentioned before, usu- 
ally it is assumed that each individual object in such an arrangement has 
a simple description - for instance, they are semi-algebraic sets defined by 
few polynomials of fixed degree. 

Arrangements of quadratic surfaces, or quadrics, in three dimensional 
space are of particular interest since they are widely used in CAD/CAM and 
computer graphics as well as in robotics ( |68j ) and computational physics 
( |58L I64j ). Therefore, it is often necessary to compute or characterize the 
intersection of quadratic surfaces and many approaches have already been 
proposed (see [MllSailllEOlESlSOllalllSTlE^ In particular, comput- 
ing the real intersection of three quadrics is an important subject in com- 
putational geometry and computer-aided geometric design (for instance, see 
[32^ [84l[83l[69]^. 
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Chionh, Goldman and Miller [32j used Macaulay's multivariate resultant 
to solve the problem in the case of finitely many intersection points. But, as 
pointed out by Xu, Wang, Chen and Sun [84], one can produce quite gen- 
eral examples where the real intersection cannot be computed using this ap- 
proach. In |84j , the computation of the real intersection of three quadrics is 
reduced to computing the real intersection of two planar curves obtained by 
Levin's method. Though useful for curve tracing. Levin's method ( |56| , 157] ) 
and its improvement by Wang, Goldman and Tu |80] has serious limitations. 
First of all, it produces a parameterization of the real intersection curve of 
two quadrics with a square-root function but does not yield information 
about reducibility or singularity of the real intersection. Secondly, Levin's 
method and similar methods ( |38^ 155] ) for computing parameterization for 
the intersection set are restricted to quadratic surfaces since higher degree 
intersection curves cannot be parameterized easily. 

In another direction, Chazelle, Edelsbrunner, Guibas and Sharir [31] 
showed how to decompose an arrangement of m objects in R*^ into 0*(m^*^~^) 
simple pieces. This was further improved by Koltun in the case k = 4 [54] . 
However, these decompositions while suitable for many applications, are not 
useful for computing topological properties of the arrangements, since they 
fail to produce a cell complex. Furthermore, arrangements of finitely many 
balls in have been studied by Edelsbrunner [39] from both combinatorial 
and topological viewpoint, motivated by applications in molecular biology. 
But these techniques use special properties of the objects, such as convexity, 
and are not applicable to general semi-algebraic sets. 

5. Review of the Results 

We review the main results of this thesis. 

5.1. Bounding the Betti Numbers. In Chapter [3] we consider the 
problem of bounding the Betti numbers, bi{S), of a semi-algebraic set S C 
R'^ defined by polynomial inequalities 

A > 0, . . . , > 0, 

where Pi G R[Xi, . . . , X^], m < k, and deg(Pj) < 2, for 1 < i < m. 
We prove (see Theorem 13. ip that for < i < A; — 1, 

min{m+l,k—i} , , , , 

, l-(^r+^. 

2 \ m J 

We first bound the Betti numbers of non-singular complete intersections of 
complex projective varieties defined by generic quadratic forms, and use this 
bound to obtain bounds in the real semi-algebraic case. Because of this new 
approach we are able to remove the constant in the exponent in the bounds 
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proved in [9|, lllj and this constitutes the main contribution which appears 
in [17] . 

5.2. Bounding the Stable Homotopy Types of a Parameterized 
Family. In Chapter H] we consider the following problem. Let 

V = {Pi,...,Prn}cR[Yi,...,Ye,Xi,...,Xk], 

with degy(Pi) < 2, degx{Pi) < d, 1 < i < m. Let S C R^^'' be a semi- 
algebaic set, defined by a Boolean formula without negations, whose atoms 
are of the form, P >0,P <0,P eV. Let vr : R^+'' R'' be the projection 
on the last k co-ordinates. Then the number of stable homotopy types (see 
Definition 12. 28p amongst the fibers 7r~^(x) n S is bounded by 

(2"*£fe(i)^(™'=) 

(see Theorem 14. ip . 

Our result can be seen as a follow-up to the recent work by Basu and 
Vorobjov |22| on bounding the number of homotopy types of fibers of general 
semi-algebraic maps (see Theorem 12. 42 p . However, our bound (unlike the 
one proven in |22] ) is polynomial in i for fixed m and k, which constitutes 
the main contribution and appears in |16j . Unfortunately, the exponential 
dependence on m is unavoidable (see Remark 14. 2p . 

Due to technical reasons, we only obtain a bound on the number of sta- 
ble homotopy types, rather than homotopy types. But note that the notions 
of homeomorphism type, homotopy type and stable homotopy type are each 
strictly weaker than the previous one, since two semi-algebraic sets might be 
stable homotopy equivalent, without being homotopy equivalent (see |73j . 
p. 462), and also homotopy equivalent without being homeomorphic. How- 
ever, two closed and bounded semi-algebraic sets which are stable homotopy 
equivalent have isomorphic homology groups. 

5.3. Algorithms and Their Implementations. In Chapter [5] we 
consider the problem of computing the first Betti Numbers of arrangements 
of compact objects in M.^ as well as computing the intersection of three 
quadratic surfaces in three dimensional space M^. 

5.3.1. Computing the Betti Numbers of Arrangements. In Chapter [T] we 
consider arrangements of compact objects in M'^ which are simply connected. 
This implies, in particular, that their first Betti number is zero. We describe 
an algorithm (see Algorithm 15. 2p for computing the number of connected 
components and the first Betti number of such an arrangement, along with 
its implementation. For the implementation, we restrict our attention to 
arrangements in and take for our objects the simplest possible semi- 
algebraic sets in M'^ which are topologically non-trivial - namely, each object 
is an ellipsoid defined by a single quadratic equation. Ellipsoids are simply 
connected, but with non-zero second Betti number. We also allow solid 
ellipsoids defined by a single quadratic inequality. This algorithm appears 
in [H]. 
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5.3.2. Computing the Real Intersection of Quadratic Surfaces. In Chap- 
ter [2] we consider the problem of computing the real intersection of three 
quadratic surfaces, or quadrics, defined by the quadratic polynomials Pi, 
P2 and P3 in M'^. We describe an algorithm for computing the isolated 
points and a linear graph embedded into M.^ (if the real intersection form 
a curve) representing the real intersection of the three quadrics defined by 
the three polynomials Pi, along with its prototypical implementation into 
the computer algebra system Maple (Version 9.5). For our implementation, 
we restrict our attention to quadrics with defining equation having rational 
coefficients. This algorithm appears in |52| . 



CHAPTER 2 



Mathematical Preliminaries 



1. Real Algebraic Geometry 



1.1. Some Notations. Let R be a real closed field and letC be an 
algebraic closed field containing R such that C = R[i]. For each m G N we 
will denote by [m] the set { 1 , . . . , m} . 

For X = (xi, . . . , Xfc) G R*^ and r G R, r > 0, we denote 



Bfc(x,r) = {y G R'' I ||y - x||2 < r^} (the closed ball), 
S''~^{x,r) = {y G R'' I ||y - x|p = r^} (the (k - l)-sphere). 

We omit both x and r from the notation for the unit sphere centered at the 
origin. For any polynomial P G R[Xi, . . . , A^fc], let 



where d is the total degree of P, the homogenization of P with respect 
to Xq. The polynomial P is Xj-regular if degjf.(P) = degP, i.e., if the 
polynomial P has a non-vanishing constant leading coefficient in the vari- 
able Xi. The gcd-free part of a polynomial P with respect to another 
polynomial Q is the polynomial P = P/ gcd{P, Q). A polynomial P G R[^] 
is square-free if there is no non-constant polynomial A G R[A] such that 

divides P. Equivalently, the polynomial P is square-free if and only if P 
is equal (up to a constant) to the gcd-free part of P and dP/dX. 

For any family of polynomials V = {Pi, . . . , Pm} C R[Ai, . . . , Xk], and 
S C R^, we denote by Zeic{V, S) the set of common zeros of V in S, i.e.. 



Let (j)he a Boolean formula with atoms of the form P = 0, P > 0, or P < 0, 
where P G P. We call (f) a P- formula, and the semi-algebraic set 5" C R^ 
defined by (p, a P-semi-algebraic set. 

If the Boolean formula (p contains no negations, and its atoms are of the 
form P = 0, P > 0, or P < 0, with P G P, then we call (p a P-cIosed 
formula, and the semi-algebraic set 5" C R'^ defined by (f), a P-closed 
semi-algebraic set. 




P\Xo,...,X,) = X^P{^,...,^) 



m 




i=l 
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For an element a G R introduce 








if a 


= 


sign (a) = < 


1 


if a 


> 




-1 


if a 


< 



A sign condition a on P is an element of {0, 1, —1}^- The realization of 
the sign condition a is the basic semi-algebraic set 

n{a) := {x G /\ sign(P(x)) = 

PeV 

A sign condition a is realizable if TZ{a) ^ 0. We denote by Sign('P) the set 
of realizable sign conditions on V. For a G Sign(P) we define the level of 
a as the cardinality 

#{P G P|a(P) = 0}. 
For each level p, < p < #P, we denote by Signp(P) the subset of Sign('P) 
of elements of level p. Furthermore, for a sign condition a let 

Z{a) := {x G R'^ I /\ P(x) = o}. 

Pe-p, <t(p)=o 

Finally, for any family of homogeneous polynomials Q = {Qi, . . . , Qm} C 
R[Xo, . . . we denote by Zer(Q,P^) (resp., Zer(Q,P^)) the set of com- 

mon zeros of Q in the real (resp., complex) projective space (resp., Pq) 
of dimension k. 

1.2. Infinitesimals. In Chapter[3]and[l]we will extend the ground field 
R by infinitesimal elements which are smaller than any positive element of 
R. The infinitesimals are used to deform our semi-algebraic sets such that 
we get very similar semi-algebraic sets having some additional properties. 

We denote by R(C) the real closed field of algebraic Puiseux series in ( 
with coefficients in R (see |20j for more details) . The sign of a Puiseux series 
in R(C) agrees with the sign of the coefficient of the lowest degree term in (. 
This induces a unique order on R(C) which makes C infinitesimal, i.e., ^ is 
positive and smaller than any positive element of R. Given a semi-algebraic 
set 5 in R'^, the extension of 5 to R(C)5 denoted Ext(S', R(C))) is the semi- 
algebraic subset of R(C)^ defined by the same quantifier free formula that 
defines S. The set Ext(5, R(C)) is well defined (i.e., it only depends on the 
set S and not on the quantifier free formula chosen to describe it). This 
is an easy consequence of the Tarski-Seidenberg principle (see for instance 

m)- 

We will also need the following remark about extensions which is again 
a consequence of the Tarski-Seidenberg transfer principle. 

Remark 2.1. Let S, T be two closed and bounded semi-algebraic subsets 
of R'^, and let R' be a real closed extension of R. Then S and T are semi- 
algebraically homotopy equivalent if and only if Ext(S', R') and Ext(T, R') 
are semi- algebraically homotopy equivalent. 
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1.3. Resultants and Subresultants. We recall next the notion of re- 
sultant and subresultant which will play an important role in the cylindrical 
decomposition and its applications (see Chapter ll.4p . We will define them 
and recall some of their properties which will be very helpful in our settings. 
But we will omit the details on how to compute them. We refer to |20j for 
more details on the algorithm. Nevertheless, it is worthwhile to mention 
that subresultants can be computed very efficiently in practice. 

Let K. be a field. Let P{X) and Q{X) be two polynomials in of 
positive degree p and q, p > 

P = apXP + --- + ao, Q = bgX'^ + --- + bo 

Next, we introduce the well-known Sylvester-Habicht matrix. 

Definition 2.2 (Sylvester-Habicht matrix). For < j < q, the j-th 
Sylvester-Habicht matrix of P and Q, denoted by SyHaj (P, Q), is the 

matrix whose rows are 

X'i-^-^P, ...,P,Q,..., XP~i~^Q considered as vectors in the basis 






Oo 











bo 



ao 

bo 






Under these conditions, the resultant of two polynomials P and Q is 
defined as follows. 

Definition 2.3 (Resultant). The (univariate) resultant of P and 

Q, denoted by Res(i-*, Q), is det(SyHao(P, Q))- 

The signed subresultants of P and Q will play a key role in what follows. 
For any j G {0, 1, . . . ,p}, the signed subresultant of P and Q of index j is 
the polynomial 

sResPj(P, Q) = sResjX^ + • • • + sResj^iX -|- sResj^o 

where sResj and each sResj^fc are elements of K defined as determinants 



of submatrices coming from SyHaj {P, Q) (see |20j for a precise definition) . 
Note that Res(P, Q) = sReso. 

We write sResPj(P, Q) (resp., Res(P, Q)) for the j-th subresultant (resp., 
resultant) of the polynomials P, Q £ M.[Xi, . . . , Xi^] with respect to X^. 
The j-th signed subresultant coefficient of P and Q, denoted by 
sResj(P, Q) or sResj, is the coefficient of X^ in sResPj(P, Q). 



^in the case p — q, we replace Q by apQ — bgP 



12 



2. MATHEMATICAL PRELIMINARIES 



Next, we notice that one of the main characteristics of subresultants is 
that they provide a very easy to use characterization of the greatest common 
divisor of two polynomials (see |20j for a proof). 

Theorem 2.4. Let P, Q £ R[^] be two polynomials of degree p and q. 
Then the following are equivalent: 

(1) P and Q have a gcd of degree j 

(2) sReso(P, Q) = ... = sResj_i(P, Q) = 0, sResj{P, Q) / 

In this case, sResP j{P,Q) is the greatest common divisor of P and Q. 

The following well-known theorem is very helpful. 

Theorem 2.5 (The Extension Theorem). Let P,Q £ C[Xi, Xk-i] [Xk], 

P = apiXi, Xk-i)Xl + . . . + ao(Xi, . . . , Xk-i) 

Q = 6,(Xi, . . . , Xk-i)Xl + ■■■ + boiXi, . . . , Xk-i). 

Let (xi, . . . , Xfc_i) G C*^"^ and assume that Res(P, (5)(xi, . . . , Xfc_i) = 0, 
then either 

(1) flp or bq vanish at (xi, . . . or 

(2) there is a number x^ G C such that P and Q vanish at (xi , . . . , Xyt) G 

Proof. See [35]. □ 

In other words, if we assume that Up and bq are in C, i.e., P and Q are X^- 
regular, and that P and Q do not have a common factor, then any solution 
(xi, . . . ,Xfc„i) G C'^"^ of the equation Res(P, Q) = can be extended to a 
solution (xi, . . . , Xfc) G C'^ of the polynomials P and Q. Note that we always 
can ensure that the polynomials are X^-regular by a change of coordinates, 
(see |83j for details). Moreover, the common factor can be detected a priori 
by computing the greatest common divisor of P and Q. 

The following proposition shows why resultants are very useful in our 
setting (see, for instance. Chapter 11.41 and [2]) . 

Proposition 2.6. Let Pi, P2 and P3 be three square-free and X^-regular 
polynomials in C[Xi, X2, X3] such that two of them do not have a common 
factor. Moreover, assume that the polynomials Res(Pi,P2) cind Res(Pi,P3) 
do not have a common factor, i.e. gcd(Res(Pi, P2), Res(Pi, P3)) = 1. Then 
the number of distinct roots of the system 

Pi(Xi, X2, X3) = 0, P2{Xi,X2,Xs) = 0, Ps{Xi,X2,Xs) = 

is finite. 

Proof. By [35], Chapter 3.6., Proposition 1, we know that Res(Pi,Pj) 
is in the elimination ideal (Pi, Pi)r]C[Xi, X2]. Therefore, by Proposition l2.5t 
only the solutions of the system 

(2.1) Res(Pi,P2) = Res(Pi,P3) = 
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can be extended to a solution of the equations (12. 6p . But there are only 
finitely many such solutions since gcd(Res(Pi, P2), Res(Pi, P3)) = 1. 

Hence, let (x, y) be a solution of the equations (12. ip . Then every 
Pi{x,y, X^) is not identically zero, as all of them are Xa-regular. In partic- 
ular, they only have finitely many solutions. Now, the claim follows. □ 

1.4. The Cylindrical Decomposition. 

1.4.1. Definition. One basic ingredient in most algorithms for comput- 
ing topological properties of semi-algebraic sets is an algorithm due to 
Collins |33j . called cylindrical decomposition, which decomposes a given 
semi-algebraic set into topological balls. In this chapter, we recall some 
facts about the cylindrical decomposition which can be turned into an algo- 
rithm for solving several important problems. For instance, computing the 
topology of planar curves (see Chapter 11.4.2^ , computing the (real) inter- 
section of quadratic surfaces (see Chapter [2]) , the general decision problem 
or the quantifier elimination problem (see |20j ). Moreover, cylindrical de- 
composition can be used to compute a semi-algebraic triangulation of a 
semi-algebraic set (see Chapter ll.Sp . For more details on the algorithm in 
the general case we refer to [33^ (2^ (3^ lit [20] . 

Definition 2.7. A Cylindrical Decomposition of R*^ is a sequence 
Si, . . . ,Sk, where, for each 1 < i < k, Si is a, finite partition of R* into semi- 
algebraic subsets (cells of level i), which satisfy the following properties: 

• Each cell S £ Si is either a point or an open interval. 

• For every 1 < i < k and every S £ Si there are finitely many 
continuous semi-algebraic functions 

^5,1 < • • • < ^S,ns ■■ S 

such that the cylinder 5 x R c R*''"^ (also called a stack over 
the cell S) is a disjoint union of cells of cSj+i which are: 

— either the graph of one of the functions ^sj, for j = 1, . . . , ns- 

{(x',Xj+i) G 5 X R I xj+i = C5,i(x')}, 

— or a band of the cylinder bounded from below and above by 
the graphs of the functions S,s,j and Csj+i, for j = 0, . . . ,ns, 
where we take ^5^0 = —00 and Cs,^s+i = +00. 

Note that a cylindrical decomposition has a recursive structure, i.e., the 
decomposition of R* induces a decomposition of R*^^ and vice-versa. 

Definition 2.8. Given a finite set V of polynomials in R[Xi, . . . , X^], a 
subset S of R*"' is P-invariant if every polynomial P in P has constant sign 
on S. A cylindrical decomposition of R'^ adapted to P is a cylindrical 
decomposition for which each cell in Sk is P-invariant. 

The following example illustrate the above definitions. 
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Figure 1. A cylindrical decomposition adapted to the unit 
sphere in 

Example 2.9 (Decomposition adapted to the unit sphere). Let 

5 = {(x,y,z) G R3 I x2 + y2 +z2 -1 = 0} 

(see Figure [1]). The decomposition of R (i.e., the line) consists of five cells 
of level 1 corresponding to the points —1 and 1 and the three intervals 
they define. The decomposition of R^ (i-e., the plane) consists of 13 cells 
of level 2. For instance, the two bands to the left and right of the circle, 
the two cells corresponding to the points (—1,0) and (1,0) and the cell that 
corresponds to the set = {(x, y) GR^|l<x<l,y = — Vl — x^}. 
The decomposition of R'^ consists of 25 cells of level 3. For instance, the 
two cells corresponding to the points (—1,0, 0) and (1, 0, 0) and the cell that 
corresponds to the set <S'3,2,2 = 'S'3,2 x {0}. For a more detailed description 
of this example see |20j . Shapter 5.1. 

The Cylindrical Decomposition Algorithm ( |33L I20j ) consists of two 
phases: the projection and the lifting phase. During the projection phase 
one eliminates the variables X^, . . . ,X2 by iterative use of (sub)-resultant 
computations. In the lifting phase the cells defined by these (sub)-resultants 
are used to define inductively, starting with i = 1, the cylindrical decompo- 
sition. 

One disadvantage of the Cylindrical Decomposition Algorithm is that it 
uses iterated projections (reducing the dimension by one in each step) and 
the number of polynomials (as well as the degrees) square in each step of 
the process. Thus, the complexity of performing cylindrical decomposition 
is double-exponential in the number of variables which makes it impractical 
in most cases for computing topological information. 



1. REAL ALGEBRAIC GEOMETRY 



15 



Nevertheless, we will see in the next chapters that it can be used quite 
efficiently for several important problems in low dimensions. 

1.4.2. Computing the Topology of Planer Curves. The simplest situation 
where the cylindrical decomposition method can be performed is the case of 
one single non-zero bivariate polynomial P £ R[Xi, X2] or a set of bivariate 
polynomials V C R[Xi,X2]. In particular, we are interested in the topology 
of the curve Zer(P, R^) (resp., of Zer(P,R^)), i.e., to determine a planar 
graph homeomorphic to Zer(P, R^) (resp., Zer(7^,R^)). 

We consider planar algebraic curves being in generic position which we 
define next. 




,.^''Zer(e) 



Figure 2. The polynomial P is in generic position with re- 
spect to Q 



Definition 2.10. Two square- free bi-variate polynomials Pi and P2 are 
in generic position with respect to the projection on the Xi-axis if the 
following conditions hold. 

(1) degiP,) = degx,{P^) (X2-regular), 

(2) gcd(Pi,P2) = l, 

(3) for all X G R the number of distinct (complex) roots of 

Pl(x,X2)=0, P2(X,X2)=0 

is or 1. 

In particular, a single bi-variate polynomial Pi is called in generic position 
with respect to P2 (resp., generic position) if Pi and dPi/dX2 ■ P2 are 
in generic position and, for / A G R, P2 / A • dPi/dX2 (resp., P2 = 1). 

It is worthwhile to mention that it is always possible to put a set of 
planar algebraic curves in generic position by a linear change of coordinates 
and computing the gcd-free part of each polynomial. Furthermore, two 
plane curves in generic position behalf nicely, i.e., their intersection points 
can be described using signed subresultant computations. The following 
proposition makes this precise. 
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Figure 3. The topology of Zev{P,R^) 



Figure 4. The topology of Zer(P, R^) with respect to Zer(Q,R^) 



Proposition 2.11. Let P, Q e R[Xi,X2] be two square-free polynomi- 
als in generic position. If (x, y) is an intersection point of Zer(P, R^) and 
Zer((5,R^), then there exists a unique j such that 

sReso(x) = • • • = sResj_i(x) = 0, sResj(x) / 
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1 sResjj„i(x) 
^ j sReSj (x) 

Proof. Let j be the unique integer such that sReso(x) = • • • = sResj_i(x) = 
and sResj(x) ^ 0. Then sResPj(P, Q)(x, X2) is the greatest common di- 
visor of the polynomials P(x, X2) and Q{x,X2) by Theorem 12.41 Since 
P and Q are in generic position, there is only one intersection point of P 
and Q with Xi-coordinate equal to x. In particular, y is the only root of 
sResPj(P, (5)(x, X2) and hence y = — (j • sResj(x))~'^sResjj_i(x). □ 

Gonzalez- Vega and Necula presented an algorithm TOP |45j which com- 
putes the topology of a plane curve. The TOP-algorithm takes a single 
bi-variate polynomial P as an input. While computing, it checks if the 
polynomial P is in generic position and performs a change of coordinates 
until the polynomial is in generic position. The TOP-algorithm outputs the 
topology of Zer(P, R^) as described below (see Algorithm I2.12p . 

For example, consider the curves given in Figure [2j The polynomial P 
(defining the two ellipses) is in generic position with respect to the polyno- 
mial Q (defining the dotted ellipse). The output of the TOP-algorithm is 
as in Figure [3l 

After some slight modifications one can use this algorithm for the fol- 
lowing two problems, which might occur simultaneously. 

(1) Computing the topology of a plane curve Zer(Pi,R^) with respect 
to another plane curve Zer(P2,R^), and 

(2) computing the common roots of two plane curves. 

Note that the proof presented in |45j can easily be adapted to those two 
problems, but the modified algorithm detects for the first problem whether 
or not the polynomial Pi is in generic position with respect to P2 and for 
the second one if Pi and P2 are in generic position. 

For our example considered above the modified TOP-algorithm output 
is as in Figure HI Note that 8 additional points are computed. 

Finally, we simply recall the in- and output of the TOP-algorithm which 
we later will use as a black-box in Chapter [2l and we refer the reader to 
|45|, I20j for more details. 

Algorithm 2.12 (TOP). 
Input: a square-free polynomial P G R[Xi,X2]. 
Output: the topology of the curve Zer(P, R^), described by 

• The real roots xi, . . . , x^ of Res(P, dP/dX2){Xi). We set by xq = 

— 00, Xj.+i = CO. 

• The number rrii of roots of P(x, X2) in R when x varies on (xj, Xj+i). 

• The number of roots of P(xj,X2) in R. We denote these roots 
by yi,i,---,yi,ni- 

• A number Cj < such that if (xj, Zj) is the unique critical point of 
the projection of Zer(P, C^) on the Xi-axis above Xj, Zj = yi^a- 
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1.4.3. Cell Adjacency. An important piece of information that we re- 
quire from the cyhndrical decomposition algorithm is that of ceh adjacency. 
In other words, we need to know given two cells in a set Si, whether the clo- 
sure of one intersects the other. In Example 12.91 for instance, we have that 
the cell corresponding to the point (—1,0,0) is adjacent to the cell 03^2,2- 

We need the following notation. We distinguish between the inter-stack 
cell adjacency of level i, which is the adjacency of cells of level i in two 
different stacks, and the intra-stack cell adjacency of level i, which is 
the adjacency of cells of level i within the same stack. 

Moreover, we use the following intuitive labeling of cells. 

• A cell in R, i.e., a cell in the induced decomposition (line) of the 
induced decomposition (plane), is denoted by (i), where the i ranges 
over the number of cells in the induced decomposition of R. Note 
that ii < 12 if and only if the cell [ii) "occurs to the left" of the 
cell {12). 

• A cell in R^, i.e., a cell in the induced decomposition of the plane, 
is denoted by (i, j), where i ranges over the number of cells in the 
line and the j ranges over the number of cells in the stack over the 
cell (i). Note that ji < j2 if and only if the cell (i, Ji) "occurs lower 
in the plane" than the cell (i, j2). 

• A cell in R^ is denoted by {i, j, k), where (i, j) is a cell in the induced 
decomposition of the plane and the k ranges over the number of 
cells in the stack over the cell Note that ki < k2 if and only 
if the cell {i,j,ki) "occurs lower" than the cell (i,j, ^^2). 

Furthermore, we distinguish among 0-cells, 1-cells, 2-cells and 3-cells of 
the cylindrical decomposition, that are points, graphs and cylinders bounded 
below and above by graphs. The adjacency between a ^-cell and k-cell will 
be denoted by {i, k}-adjacency. 

We illustrate the above notation on Example l2. 91 (Decomposition adapted 
to the unit sphere). 

Example 2.13 (cont.). For instance, the cell (2) and (4) correspond 
to the points —1 and 1 (in the line), whereas the cells (2,2) and (3,2) 
correspond the point ( — 1,0) and the set 

^3,2 = {(x,y) E R2 I -1< X < l,y = -TT^}. 

Moreover, the cell (2, 2, 2) corresponds to the point (—1, 0, 0) and the cell (3, 2, 2) 
corresponds to the set 5*3^2,2 = 'S'3,2 x {0}. 

While there are algorithms known for computing the cell adjacencies 
of a cylindrical decomposition of R'^ (see ((3l |4j), we will only be inter- 
ested in the cell adjacencies for a cylindrical decomposition adapted to fam- 
ily V C R[Xi,X2,X3] such that deg(P) < 2 and P is Xa-regular for every 
polynomial P £ V. 

It is worthwhile to mention that we do not need to compute all cell 
adjacencies. In our applications (see Chapter [5]) it suffices to compute the 
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{0, l}-inter-stack adjacencies which we can do by a simple combinatorial 
type approach. In other words, we determine the full adjacency information 
for the boundary of the semi-algebraic set by using the simpler structure 
induced by the quadratic polynomials which we describe next. 

Assume that the 0-cell and the 1-cell (i + l,j2) are adjacent in 

the induced decomposition of the plane. To be more precise, the 0-cell 
and the 1-cell (i + l,j2) correspond to a point and a curve segment of 
Zer(Res(Pm5 -Pii ^3)5 R^) where Pm and Pt are two input quadratic polyno- 
mials that are Xs-regular. We have the following two cases: 

Case 1: The stack over the 0-cell contains exactly one 0-cell 

{i,ji,k). Note, that the stack over 1-cell (i + 1,^2) must contain two 1- 
cells {i + l,j2,li) and {i + 1,^2)^2) (corresponding to graphs), since the 
polynomial Pm is of degree equal to 2 in the variable X^. Therefore, the 

0- cell (i, ji, k) must be adjacent to both cells (i + 1, j2, ^i) and {i + 1, j2, ^2)1 
since the semi-algebraic set Si is closed. 

Case 2: The stack over the 0-cell (i,ji) contains two 0-cells {i,ji,ki) 
and ^2)- As above, the stack over the 1-cell (2 + 1,^2) must contain two 

1- cells {i+l,j2,h) and (i+1, j2, ^2)- Remember that both stacks are ordered 
from the bottom to the top. Hence, the cells fci) and (i + 1,^2, ^1) as 
well as the cells {i,ji,k2) and {i + 1,^2,^2) must be adjacent for the same 
reason as above. It is worthwhile to mention that is not possible to have just 
one 1-cell above (i + 1, j2)) i-e., ii = i2, by the properties of the cylindrical 
decomposition. 

1.5. Triangulation of Semi-algebraic Sets. Another important prop- 
erty of closed and bounded semi-algebraic sets is that they are homeomor- 
phic to a simplicial complex. The following makes this statement precise. 

Let ao,...,ap be points of R'^ that are affinely independent. The p- 
simplex with vertices ao, . . . ,ap is 

p 

[ao, ...,ap] = {Aoao H h XpUp \ ^ Aj = 1 and Ao, . . . , Ap > 0} 

i=0 

Note that the dimension of [ao, . . . , ap] is p. 

An q-face of the p-simplex s = [oq ap] is any simplex s' = [bo, ■ ■ ■ , bg] 
such that 

{bo, ...,bq}c {ao, . . . ,ap} 

The open simplex, denoted by s", corresponding to a simplex s consists of 
all points of s which do not belong to any proper face of s: 

p 

s° = {ao, ... ,ap) = {Aoao H h Apap | ^ Aj = 1 and Ao > 0, . . . , Ap > 0} 

i=0 

A simplicial complex K in R'^ is a finite set of simplices in R'^ such that 
s, s' G K implies 

• every face of s is in K, 
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• s n s' is a common face of both s and s'. 
A triangulation of a semi-algebraic set S" is a simplicial complex K together 
with a semi- algebraic homeomorphism /i : |K| — > S, where the set |K| = 
UseK ^ realization of K. 

A triangulation of S respecting a finite family of semi-algebraic 
sets Si, . . . , Sn contained in 5" is a triangulation (K, h) such that each Sj is 
the union of images by h of open simplices of K. 

We have the following theorem. 

Theorem 2.14. Let S CK^ be a closed and bounded semi- algebraic set, 
and let Si, . . . ,Sn be semi- algebraic subsets of S. There exists a triangulation 
of S respecting Si, . . . ,Sn- Moreover, the vertices of K can be chosen with 
rational coefficients. 

Proof. See [20] □ 

For example, let S be a closed and bounded subset of R'^ such that 
S = ljr=i C R^. Then Theorem 12.141 implies that there is a triangulation 
(K, h) of S such that for every simplex s E K and 1 < i < n either h{s)nSi = 
h{s) or h{s) nSi = $. 

Finally, note that one can compute a triangulation of a closed and 
bounded semi-algebraic set using the cylindrical decomposition which de- 
composes a given semi-algebraic set into double exponential number (in the 
dimension) of topological balls. 

1.6. Triviality of Semi-algebraic Mappings. The finiteness of the 
topological types of algebraic subsets of R'^ defined by polynomials of fixed 
degree is an easy consequence of Hardt's triviality theorem, which we recall 
next. 

Theorem 2.15 (Hardt's triviality theorem [491, HQ]). Let 5 C R" and 

T C R^ 6e semi- algebraic sets. Given a continues semi- algebraic function f : 
S ^ T, there exists a finite partition of T into semi- algebraic sets T = 
Uig/ ^0 that for each i and any G Ti, Tj x /~^(xj) is semi- algebraically 
homeomorphic to f~^{Ti). 

Hardt's theorem is a corollary of the existence of cylindrical decompo- 
sitions (see Chapter II. 4p . which implies a double exponential (in n) upper 
bound on the cardinality of the set I. Moreover, it follows that one can 
always retract a closed semi-algebraic set to a closed and bounded set. The 
following proposition makes this precise. 

Proposition 2.16 (Conic structure at infinity). Let S cR'' be a closed 
semi- algebraic set. There exists r € H, r > 0, such that for every r' , r' > r, 
there is a semi- algebraic deformation retraction from S to Sr' = 5nBfc(0, r') 
and a semi- algebraic deformation retraction from Sj.' to Sr. 



Proof. See [20], Proposition 5.49. 
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2. Algebraic Topology 

2.1. Some Notations. In this chapter we recall the basic objects from 
algebraic topology like homology and co-homology theory. Unless otherwise 
noted, we will consider vector spaces over Q in what follows next. 

Given a simplicial complex K, we denote by Cp(K) the vector space 
generated by the p-dimensional oriented simplices of K. The elements of 
Cp(K) are called the p-chains of K. For p < 0, we define Cp(K) = 0. 

Given an oriented p-simplex s = [oq, . . . , Op], p > 0, the boundary of s is 
the {p — l)-chain 

9p{s) = ^ {-iy[ao,- ■ ■ ,ai_i,ai,ai+i,. . . ,ap\, 

0<i<p 

where fij means that the is omitted. For p < 0, we define dp = 0. The 
map dp extends linearly to a homomorphism 

dp : Cp(K) ^ Cp_i(K). 

Thus, we have the following sequence of vector space homomorphism with 

dp-i odp = 0, 

> Cp{K) ^ Cp„i(K) ^ Cp_2(K) ^ . . . A, Co(K) ^ 

The sequence of pairs {(Cp(K), 9p)}pgN, denoted by C,(K), is called the 
simplicial chain complex. 

We denote by Hp(K) the p-th simplicial homology group of K, that 

is 

Hp(C.(K)) = Zp(C.(K))/Bp(C.(K)), 

where Zp(C,(K)) = Ker((9p) is the subspace of p-cycles, and Bp(C,(K)) = 
Im(5p4-i) is the subspace of p-boundaries. 

Note that Hp(K) is a finite dimensional vector space. The dimension of 
Hp(K) as a vector space is called the p-th Betti number of K and denoted 
by 6p(K). We will denote by b{K) the sum Ylp>o^p(^)- 

Next, we define the dual notion of cohomology groups. 
We denote by C{K) = Hom(Cp(i^), Q) the vector space dual to Cp{K), and 
by 6^ the co-boundary map 6P : CP{K) CP+^K) which is the homomor- 
phism dual to 5p+i in the simplicial chain complex C,{K). More precisely, 
given LO G C^(K), and a p + 1-simplex [oq, . . . , flp+i] of K, then 

(5w([ao, . . . ,ap+i]) = ^ (-l)'u;([ao,...,ai_i,dj,ai+i,...,ap+i]) 

0<j<p+l 

Thus, we have the following sequence of (dual) vector space homomorphism, 

^ C°(K) ^ &{K) C2(K) ^ . . . CP(K) ^ CP+\K) ^ • • • , 

with (5^+^0(5^ = 0. The sequence of pairs {(C''^(K), 5P)}pgN, denoted by 
C*{K), is called the simplicial cochain complex. 
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We denote by W{K) the p-th simplicial cohomology group of K, 

that is 



HP(C'(K)) = ZP(C'(K))/BP(C*(K)), 
where ZP(C'(K)) = Ker(9P-i) is the subspace of p-cocycles, and BP(C*(K)) = 



lm{dp) is the subspace of p-coboundaries. 

Note that H^(K) is a finite dimensional vector space and its dimension as 
a vector space is equal to bp(K) . To be more precise, we have by the Universal 
Coefficient Theorem for cohomology (see |51| . Theorem 3.2, page 195) that 
H^(C*(K)) and Hp(C,(K)) are isomorphic for every p > 0. Moreover, the 
cohomology group H^(K) can be identified with the vector space of locally 
constant functions on |K| (see |20j . Proposition 6.5). 

Next, we define simplicial (co)-homology groups for a closed semi- algebraic 
set. Let S* C R'^ be a closed semi-algebraic set. By Proposition 12.161 (Conic 
structure at infinity), there exists r S R, r > 0, such that for every r', 
r' > r, there is a semi-algebraic deformation from S to S^' =30 Bfc(0, r') 
and a semi- algebraic deformation from Sr' to Sr- Note that the set Sr is 
closed and bounded. By Theorem 12.141 the set Sr can be triangulated by 
a simplicial complex K with rational coordinates. Choose a semi-algebraic 
triangulation / : |K| — > 5^, then for p > the homology groups Hp(S) 
are Hp(K) (resp., cohomology groups HP(S) are H^(K)). Note that the 
(co)-homology groups do not depend on the particular triangulation. The 
dimension of Hp(S) as a vector space is called the p-th Betti number of 
S and denoted by bp{S). We will denote by b{S) the sum X]p>o bp{S). 

For completeness we now consider a basic locally closed semi-algebraic 
set S which is, by definition, the intersection of a closed semi-algebraic set 
with a basic open one. Let S be the (one point) Alexandroff compactification 
of S. Then the dimension of 11^(5") as a vector space is called the p-th Betti 
number of S and denoted by bp{S). This definition is well-defined since the 
Alexandroff compactification 5 of 5* is closed, bounded, unique (up to semi- 
algebraic homeomorphism) and semi-algebraically homeomorphic to S. We 
will denote by b{S) the sum X]p>o ^p('S'). Note that the homology groups 
of a semi-algebraic set S C R^ are finitely generated. Hence, the Betti 
numbers bi(S) are finite. 

We illustrate Betti numbers with the following example. 




Figure 5. The hollow torus 



2. ALGEBRAIC TOPOLOGY 



23 



Example 2.17. Let S be the hollow torus in R'^ (see Figure [5]), then 
6o(5) = l, bi{S) = 2, b2iS) = 1 and bp{S) =0,p> 2. 

Intuitively, bp{S) measures the number of p-dimensional holes in the set 
S. The zero-th Betti number, bo^S), is the number of connected components. 

Similarly, one can define 6^(5, Z2), the p-th Betti number with coeffi- 
cients in Z2, as the Z2-vector space dimension of Hp(S', Z2). We denote by 
6(5,^2) the sum Ep>o ^p('S'i ^2)- 

It follows from the Universal Coefficients Theorem, that 

k{S,Z2) > bi{S) 

(see [51], Corollary 3.A6 (b)). 

Hence, any bounds proved for Betti numbers with Z2-coefficients also 
apply to the ordinary Betti numbers (with coefficients in Q). 

2.2. The Mayer- Vietoris Theorem. We have seen in Chapter 11.41 
that we can use the cylindrical decomposition in order to decompose a 
semi- algebraic set into smaller pieces. The Mayer- Vietoris inequalities (see 
Proposition 12.19]) bound the Betti numbers of the union (resp., intersection) 
of semi-algebraic sets in terms of intersections (resp., unions) of fewer semi- 
algebraic sets. This will be very useful in Chapter [3] and Chapter [H We 
first recall a semi-algebraic version of the Mayer- Vietoris theorem. 

Theorem 2.18 (Semi-algebraic Mayer- Victor is). Let Si and S2 be two 
closed and bounded semi- algebraic subsets ofK^. Then there is a long exact 
sequence. 

> Hp(5inS2) ^ Hp(5i)eHp(52) ^ Hp(SiU52) ^ Hp_i(5in52) ^ • • • 

Proof. By Theorem 12.141 there is a triangulation of Si U 5*2 that is 
simultaneously a triangulation of ^i, ^2, and S'in52. Let Ki be the simplicial 
complex corresponding to Si. Then there is a a short exact sequence of 
simplicial chain complexes, 

^ c.{Ki n K2) ^ c.(Ki) e c,{K2) c,{Ki u ^2) ^ o 

The claim follows by a standard argument about short and long exact se- 
quences (see |20j . Lemma 6.10). □ 

From the exactness of the Mayer- Vieotoris sequence, we have the follow- 
ing proposition. 

Proposition 2.19 (Mayer- Vietoris inequalities). Let be Si, . . . , Sn sub- 
sets of R'^ be all open or all closed. Then for each i > we have, 

(2-2) hi u sA<Y^ I n 

\l<j<n I Jc[n] \jeJ 
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and 

(2-3) bj n SA<Y1 biM#J)~i 

\l<j<n / JC[n] 

Proof. Follows from [20], Proposition 7.33. □ 

The following proposition characterizes 6o cind bi in a special case of 
unions of simplicial complexes. It is a slightly strengthened version of a 
similar proposition appearing in [211, 120j . We do not require that the com- 
plexes Ai be acyclic, but only that their first co-homology group vanishes. 
We need the following notations. 

Let Ai, . . . , An be sub-complexes of a finite simplicial complex A such 
that 

• each Ai is connected, i.e., H°(^i) = Q, 

• A = Ur=i -4i, and 

• R^{Ai) = 0, I <i <n. 

Note that the intersections of any number of the sub-complexes, Ai, is again 
a sub-complex of A. We will denote by Aij the sub-complex Ai n Aj, and 
by Aij^i the sub-complex Ai n Aj n Ai. 

Recall that H^(K) can be identified as the vector space of locally constant 
functions on the simplicial complex K. Hence, we can define the following 
sequence of generalized restriction homomorphisms. 

Let (/.eei<,<,HO(A,), define 

{So(p)i,j = </>iUi,, - 4>j\A,j 
and let V G 0i<i<j<„, H°(^ij), define 

We now are able to state our proposition. 

Proposition 2.20. Let Ai, . . . , An be sub- complexes of a finite simplicial 
complex A such that A = IJiLi ^« '^'^'^ Z'"" each i,l < i < n, 

(1) R^{Ai) = Q, and 

(2) R\A,) = 0. 

Let the homomorphisms 6q and Si in the following sequence be defined as 
above. 

i i<j i<j<i 

Then, 

(1) bo{A) = dim(Ker((^o)), 

(2) bi{A) = dim(Ker((^i)) - dim(Im((5o)). 

Proof. Follows from |20], Theorem 6.9. □ 
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Remark 2.21. One could use the so-cahed generahzed Mayer- Vietoris 
sequence and some spectral sequence argument in order to prove Proposi- 
tion I2.20[ We refer to \10\ I15j for more details. 

2.3. Smith Theory. Li Chapter[3]we will reduce the problem of bound- 
ing the Betti numbers of a semi-algebraic set to the problem of bounding 
the Betti numbers of some real projective algebraic sets. Using the Smith 
inequality (see Theorem 12. 221 below) allows us to relate the Betti numbers of 
these real projective algebraic sets to the corresponding complex projective 
algebraic sets. As we will see in Chapter 13.21 we have precise information 
about the corresponding complex projective algebraic set. Before we recall 
a version of the Smith inequality, we need the following. 

Let X be a compact topological space and c : X ^ X an involution. We 
regard X as a G-space, where G = {id,c} = Z2. We denote by X' = X/c 
the orbit space, and hy F = Fix c, the fixed point set of the involution c. 
Moreover, we identify F with its image in X' . 

Then there are two exact sequences, called (homology and cohomol- 
ogy) Smith sequences of (X,c): 

> Hp+i(X',F;Z2) ^ Hp(X',F;Z2)eHp(F;Z2) ^ Rp{X;Z2) ^ Rp{X',F;Z2) 

.Hp(X',F;Z2) ^HP(X;Z2) ^Hf(X',F;Z2)eHP(F;Z2) ^HP+i(X',F;Z2) 

We refer the reader to |78j . p. 131, for more details. 

Next, we state a version of the Smith inequality which follows from the 
exactness of the Smith sequence. We consider the special case where X is a 
complex projective algebraic set defined by real forms, with the involution 
taken to be complex conjugation. Then we have the following theorem. 

Theorem 2.22 (Smith inequality). Let Q C R[Xo, . . . ,Xk] be a family 
of homogeneous polynomials. Then, 

6(Zer(Q,P|),Z2) < 6(Zer(Q, P^), Z2). 

2.4. Alexander Duality. In Chapter [3l we also use the well-known 
Alexander duality theorem which relates the Betti numbers of a compact 
subset of a sphere to those of its complement. 

Theorem 2.23 (Alexander Duality). Let r > 0. For any closed sub- 
set A C S''{0,r), 

RiiS''{0,r)\A)^R''-'-HA), 
where W{A), < i < k — 1, denotes the reduced cohomology group of A. 

Proof. See [62], Theorem 6.6. □ 

2.5. The Betti Numbers of a Double Cover. Let X be a topolog- 
ical space. A covering space of X is a space X together with a continuous 
surjective map f : X ^ X, such that for every x £ X there exists an open 
neighborhood U oi x such that f~^{U) is a disjoint union of open sets in X 
each of which is mapped homeomorphically onto U by /. In particular, if 
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for every x E X the fiber /~^(x) has two elements, we speak of a double 
cover. 

The following proposition relates the Betti numbers (with Z2 coefficients) 
of a finite simplicial complex to its double cover. Note that the proposition 
is no longer true for Betti numbers (with Q-coefficients). A simple coun- 
terexample is provided by the 2-torus which is a double cover of the Klein 
bottle, for which the stated inequality is not true for i = 2 for Betti numbers 
(with Q-coefficients) . 

Proposition 2.24. Let X be a finite simplicial complex and f : X ^ X 
a double cover of X. Then for each i > 0, 

h{X,Z2) <2bi{X,Z2). 

Proof. Let 

:C.(X,Z2) ^C.(X,Z2) 

denote the chain map sending each simplex of X to the sum of its two 
preimages in X. Let 

V'. :C.(X,Z2) ^C.(X,Z2) 

be the chain map induced by the covering map /. 

It is an easy exercise to check that the following sequence is exact, 

C.(X,Z2) ^ C.(X,Z2) ^ C.(X,Z2) 0. 

The corresponding long exact sequence in homology, 

> Ri{X,Z2) ^ H,(X,Z2) ^ Hi(X,Z2) • • • 

gives the required inequality. □ 

Remark 2.25. The above proof is due to Michel Coste. 

2.6. The Betti Numbers of a Projection. The following proposi- 
tion gives a bound on the Betti numbers of the projection n{S) of a closed 
and bounded semi-algebraic set S in terms of the number and degrees of 
polynomials defining S. 

Proposition 2.26 ([43]). LetR be a real closed field and letir : R™+*^ 
be the projection map on to last k co-ordinates. Let S C R™+'^ be a closed 
and bounded semi- algebraic set defined by a Boolean formula with s distinct 
polynomials of degrees not exceeding d. Then the n-th Betti number of the 
projection 

b„(^(5)) < (mn(i)0('=+"'"). 
Proof. See HS. □ 



2.7. The Smale-Vietoris Theorem. In Chapter H] we also need the 
following version of the well-known Smale-Vietoris Theorem |71j . 

Theorem 2.27 ( |71j ). Let S andT be closed and bounded semi-algebraic 
sets, and f : S ^ T a continuous semi- algebraic map such that f~^{y) is 
contractible for every y G T. Then the map f is a homotopy equivalence. 
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2.8. Stable homotopy equivalence and Spanier- Whitehead du- 
ality. For any finite CW-complex X we will denote by S{X) the suspension 
of X, which is the quotient of X x [0, 1] by collapsing X x {0} to one point 
and X X {1} to another point. 

Recall from [72] that for two finite CW-complexes X and Y, an element 

of 

(2.4) {X;y} = lim[S^(X),S*(y)] 

i 

is called an S-map (or map in the suspension category). (When the context 
is clear we will sometime denote an S-map / G {X; Y} by / : X — > y). 

Definition 2.28. An S-map / G {X;Y} is an S-equivalence (also 
called a stable homotopy equivalence) if it admits an inverse G 
In this case we say that X and Y are stable homotopy equiva- 
lent. 

If / G {X;y} is an S-map, then / induces a homomorphism, 

A :H,(X) ^H,(y). 

The following theorem characterizes stable homotopy equivalence in terms 
of homology. 

Theorem 2.29. [73] Let X and Y be two finite CW-complexes. Then 
X and Y are stable homotopy equivalent if and only if there exists an S-map 
f G {X;y} which induces isomorphisms /* : Ilj(X) — > Hj(y) (see [36] . pp. 
604) for all i > 0. 

In order to compare the complements of closed and bounded semi- 
algebraic sets which are homotopy equivalent, we will use the duality the- 
ory due to Spanier and Whitehead [72] . We will need the following facts 
about Spanier- Whitehead duality (see [36] . pp. 603 for more details). Let 
X C S" be a finite CW-complex. Then there exists a dual complex, denoted 
DnX C S"'\X. The dual complex DnX is defined only upto S-equivalence. 
In particular, any deformation retract of S"'\X represents DnX. Moreover, 
the functor D„ has the following property. If y C S" is another finite CW- 
complex, and the S-map represented hy <j) : X ^ Y is a stable homotopy 
equivalence, then there exists a stable homotopy equivalence -Dn^- More- 
over, if the map (p : X ^ Y is an inclusion, then the dual S-map Dn(j) is also 
represented by a corresponding inclusion. 

Remark 2.30. Note that, since Spanier- Whitehead duality theory deals 
only with finite polyhedra over M, it extends without difficulty to general 
real closed fields using the Tarski-Seidenberg transfer principle. 

2.9. Homotopy colimits. Let A = {Ai, . . . ,An}, where each Ai is a 
sub-complex of a finite CW-complex. 
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Let A[„] denote the standard simplex of dimension n — 1 with vertices 
in [n]. For I C [n], we denote by A/ the (#/ — l)-dimensional face of A^^j 

corresponding to /, and by Aj the CW-complex ^ Ai. 

The homotopy cohmit, hocohm(^), is a CW-complex defined as follows. 
Definition 2.31. 

hocolim(^) = [J A/ X A// ~ 

/C[n] 

where the equivalence relation ~ is defined as follows. 

For / C J C [n], let s/.j : A/ ^ Aj denote the inclusion map of the 
face A/ in Aj, and let ii j : Aj ^ Ai denote the inclusion map of Aj in 
Ai. 

Given (s,x) G A/ x A/ and (t,y) G Aj x Aj with I C J, then (s,x) ~ 
(t,y) if and only if t = s/,j(s) and x = i/,j(y). 

We have a obvious map 

/ : hocolim(^) — > colim(^) = A, 

ie[n] 

sending (s, x) i— > x. It is a consequence of the Smale-Vietoris theorem (see 
Theorem ??) that 

Lemma 2.32. The map 

f : hocolim(^) — > colim(^) = 

i&[n] 

is a homotopy equivalence. 

Now let A = {Ai, . . . , An} (resp. B = {Bi, . . . , Bn}) be a set of sub- 
complexes of a finite CW-complex. For each / C [n] let // G {Ai;Bj} be a 
stable homotopy equivalence, having the property that for each / C J C [n], 
fj = fi\Aj- Then we have an induced S-map, / G {hocolim(^); hocolim(S)}, 
and we have that 

Lemma 2.33. The induced S-map f G {hocolim(^); hocolim(;S)} is a 
stable homotopy equivalence. 

Proof. Using the Mayer- Vietoris exact sequence it is easy to see that 
if the //'s induce isomorphisms in homology, so does the map /. Now apply 
Theorem [2291 □ 

3. The Topology of Algebraic and Semi-Algbraic Sets 

3.1. Bounds on the Topology of Semi- Algebraic Sets. The initial 
result on bounding the Betti numbers of semi- algebraic sets defined by poly- 
nomial inequalities was proved independently by Oleinik and Petrovskii |65j , 
Thom [76] and Milnor [63]. They proved: 
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Theorem 2.34. [65l EHl [63] Let 

7' = {Pi,...,P„} cR[Xi,...,Xfc] 
with deg(Pi) < d, 1 < i < m and let S C R'^ he the set defined by 

A > 0, . . . , > 0. 

Then 

h{S) = 0{mdf. 

Notice that the theorem includes the case where the set S is a real 
algebraic set. Moreover, the above bound is exponential in k and this ex- 
ponential dependence is unavoidable (see Example 12.351 below) . Recently, 
the above bound was extended to more general classes of semi- algebraic sets. 
For example, Basu [llj improved the bound of the individual Betti numbers 
of P-closed semi-algebraic sets while Gabrielov and Vorobjov |44j extended 
the above bound to any P-semi-algebraic set. They proved a bound of 
0{m'^df. Mor cover, Basu, Pollack and Roy [19] proved a similar bound for 
the individual Betti numbers of the realizations of sign conditions. 

Example 2.35. The set S (ZR^ defined by 

Xi(^i-l) >0,...,Xfc(Xfc-l) >0, 

has 6o(S) = 2^ 

However, it turns out that for a semi-algebraic set S d R^ defined by 
m quadratic inequalities, it is possible to obtain upper bounds on the Betti 
numbers of S which are polynomial in k and exponential only in m. The 
first such result was proved by Barvinok who proved the following theorem. 

Theorem 2.36. ^ Let S d R^ he defined hy 

Pi > 0, . . . , P„ > 0, 

with deg(Pi) < 2, 1 < i < m. Then, h{S) < k^'-'^\ 

Theorem 12.361 is proved using a duality argument that interchanges the 
roles of k and m, and reduces the original problem to that of bounding the 
Betti numbers of a semi-algebraic set in R* defined by k^^^^ polynomials 
of degree at most k. One can then use Theorem 12.341 to obtain a bound of 
j^O{m) ^ rpj^g constant hidden in the exponent of the above bound is at least 
two. Also, the bound in Theorem l2.36l is polynomial in k but exponential in 
m. The exponential dependence on m is unavoidable as remarked in [9], but 
the implied constant (which is at least two) in the exponent of Barvinok's 
bound is not optimal. 

Using Barvinok's result, as well as inequalities derived from the Mayer- 
Vietoris sequence, Basu proved a polynomial bound (polynomial both in k 
and m) on the top few Betti numbers of a set defined by quadratic inequal- 
ities. More precisely, he proved the following theorem. 



30 



2. MATHEMATICAL PRELIMINARIES 



Theorem 2.37. [IT] Let £ > and let S C R'' be defined by 
with deg(Pj) < 2. Then 

Notice that for fixed i, the bound in Theorem l2.37l is polynomial in both 
m and k. 

3.2. Bounds on the Topology of Complex Algebraic Sets. By 

separating the real and imaginary parts one can consider a complex algebraic 
set X C as a real algebraic subset of R^'^. Unfortunately, real and 
complex algebraic sets do not have the same properties. To be more precise, 
an irreducible algebraic subset of C*^ having complex dimension n, considered 
as an algebraic subset of R^'^ is connected, not bounded (unless it is a point) 
and has local real dimension 2n at every point (see, for instance, |27] ). But 
this is no longer true for real algebraic sets as we will see in the following 
examples. 

Example 2.38 ([22]). (1) The circle {(x,y) G R^ | x^ + y2 = 1} 
is bounded. 

(2) The cubic curve {(x, y) G R^ | x^ + y^ — x'^ = 0} has an isolated 
point at the origin. 

However, in Chapter[3]we will show how to reduce the problem of bound- 
ing the Betti numbers of a real algebraic set to the problem of bounding 
the Betti numbers of a complex projective algebraic set involving the same 
polynomials. Moreover, this complex projective algebraic set will have the 
property that is a non-singular complete intersection, which we define next. 

Definition 2.39. A projective algebraic set X C of codimension n 
is a non-singular complete intersection if it is the intersection of n non- 
singular hypersurfaces in Fq that meet transversally at each point of the 
intersection. 

Next, we recall some results about the Betti numbers of a complex pro- 
jective algebraic set which is a non-singular complete intersection. We need 
the following notation. 

Fix a j-tuple of natural numbers d = {di, . . . ,dj). Let Xq = Zer({Qi, . . . , Qj}, P^), 
such that the degree of Qi is di, denote a complex projective algebraic set 
of codimension j which is a non-singular complete intersection. 

Let b{j, k, d) denote the sum of the Betti numbers with Z2 coefficients of 
Xq. This is well defined since the Betti numbers only depend only on the 
degree sequence and not on the specific Ac (see, for instance, |41j ). 

The function b{j,k,d) satisfies the following (see |26j ): 
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b{j, k, d) 

where 



c{j, k, d) a k — j is even, 

2{k — j + 1) — c{j, k, d) M k — j is odd, 



A; + 1 if J = 0, 

c(j, k,d) = ^di... dj if j = k, 

dkc{j -l,k-l, (di,.. . ,4-i)) - (4 - l)c(i,/c - l,d) ifO < j <k. 

In the special case when each 4 = 2, we denote by b{j, k) = b{j, k,(2, . . . , 2)). 
We then have the following recurrence for b{j,k). 



b{j,k) 

where 



q{j, k) \i k — j is even, 

2{k — i + 1) — q{j^ k) \i k — j is odd, 



Next, we show some properties of q{j, k). 



k + 1 if J = 0, 

2^ iij = k, 

[2q{j - 1, A; - 1) - q{j, k - 1) if < j < k. 



Lemma 2.40. (1) q{l,k) = k+l/2{l- {-!)'') andq{2,k) = {-l)''k+ 
k. 

(2) For2<j<k, |g(j,A;)| <2^-i(/J. 

(3) For2<j <k and k-j odd, 2{k - j + 1) - q{j, k) < 2^'^ J . 

Proof. The first part is shown by two easy computations and noting 
that 

2(A: - 2 + 1) - q{2, k) = 2k - 2 if k - 2 is odd. 

Hence, we can assume that the statements are true for k — 1 and that 3 < j < 
k. Note that for the special case j = /c — 1, we have that 2^"""^ < 2*''~^(^~2) 
since k > 2. Then 

\qU,k)\ = \2q{j -l,k-l)-q{j,k-l)\ 
< 2\qij-l,k-l)\ + \q{j,k-l)\ 



2^ 



~i 



k 



32 



2. MATHEMATICAL PRELIMINARIES 



and, for k — j odd, 

2(k-j + l)-q{j,k) = 2{k-j + l)-2q{j-l,k-l)+q{j,k-l) 

< \2i{k - 1) - (j - 1) + 1) - q{j - 1, A; - 1)1 
+ \q{j-l,k-l)\ + \qij,k-l)\ 



j-2j \j-lJJ 

□ 

Hence, we get the following bound for 6(j, k). 
Theorem 2.41. 

(1) 6(1 A;) = i ^''^'^ ~ «/A:ise?;en, 

' I g(0, k) if k is odd, 

(2) bij,k)<2^-^{.'_^),for2<j<k. 

Proof. Follows from Lemma [2ll0l □ 

3.3. Bounds on the Topology of Parametrized Semi-algebraic 
Sets. Let vr : R^"'"'^ — > R'^ be the projection map on the last k co-ordinates, 
and for any S C R^^'^ we will denote by its the restriction of vr to S. 
Moreover, when the map vr is clear from context, for any x G R'^ we will 
denote by the fiber vr~^(x) n S. One way to interpret this setting is 
that the set S depends on k parameters and vr is the projection onto the 
parameter space. 

Hardt's triviality theorem (see Theorem 12.15^ implies that there ex- 
ists a semi-algebraic partition {Tjjjg/ of R'^ having the following property. 
For each i £ I and any point x G Tj, the pre-image vr~^(rj) n 5" is semi- 
algebraically homeomorphic to x Tj by a fiber preserving homeomorphism. 
In particular, for each i £ I, all fibers (Sx, x G Tj are semi-algebraically 
homeomorphic. 

As mentioned in Chapter [L6] the existence of cylindrical decompositions 
implies a double exponential (in k and i) upper bound on the cardinality of 
/ and, hence, on the number of homeomorphism types of the fibers of the 
map TTs- No better bounds than the double exponential bound are known, 
even though it seems reasonable to conjecture a single exponential upper 
bound on the number of homeomorphism types of the fibers of the map vr^. 

In |22j . Basu and Vorobjov considered the weaker problem of bounding 
the number of distinct homotopy types occurring amongst the set of all 
fibers of Sx, and they proved a single exponential upper bound (in k and £) 
on the number of homotopy types of such fibers. 

They proved the following theorem. 
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Theorem 2.42. ^ LetV C R[Yi, . . . ,Yi,Xi, . . . ,Xk], with deg(P) < d 
for each P £ V and cardinality = m. Then there exists a finite set 
^4 C R'' with 

#A < (2^mfcd)«(^^^) 

such that for every x G R'^, there exists z €z A such that for every V-semi- 
algebraic set S C R^"*"^, the set is semi-algebraically homotopy equivalent 
to Sz- In particular, for any fixed V -semi- algebraic set S, the number of 
different homotopy types of fibers 5'x for various x G vr(S') is also bounded 
by 

{2'mkdf^^'\ 

Notice that the bound in Theorem 12.421 is single exponential in kL The 
following example, which also appears in |22j, shows that the single expo- 
nential dependence on £ is unavoidable. 

Example 2.43. Let P G RfFi, . . . ,Y(\ ^ K\Yi, . . . , y^, X] be the poly- 
nomial defined by 

I d-i 

i=l j=0 

The algebraic set defined by P = in R^"*"^ with coordinates Yi, . . . ,Ye, X , 
consists of d^ lines all parallel to the X axis. Consider now the semi-algebraic 
set S C R^+^ defined by 

(P = 0) A {0 < X <Yi + dY2 + d^Ys + ■ ■ ■ + d^-^Y^). 

It is easy to verify that, if tt : R^"^^ — > R is the projection map on the 
X co-ordinate, then the fibers 5x, for x G {0, 1, 2, . . . , — 1} C R are 0- 
dimensional and of different cardinality, and hence have different homotopy 
types. 

3.4. Some Useful Constructions. In this chapter, we recall some 
very useful constructions for semi-algebraic subsets of R*^ which are well- 
known in real algebraic geometry. 

Let = {Pi, . . . , P„} C R[Xi, . . . , Xfc] with deg(Pi) < 2, 1 < i < m. 
Let 5 C R'' be the basic semi-algebraic set defined by 

5 = {x G R'^ I Pi (x) > 0, . . . , P^(x) > 0}. 

Let 1 ^ e > be an infinitesimal, and let 

fc 

Pm+l = 1 ~ ^i- 
i=l 

Let ^6 C R(e) be the basic semi-algebraic set defined by 

= {x G R(e)'^ I Pi(x) > 0, . . . , P^(x) > 0, P™+i(x) > 0}. 
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2. MATHEMATICAL PRELIMINARIES 



Proposition 2.44. The bounded set Sh and the set Ext(5, R(e)) are 
homotopy equivalent. Moreover, the homology groups of the Sb and S are 
isomorphic. 

Proof. It follows from Proposition 12.161 (Conic structure at infinity) 
that the semi-algebraic set Si, has the same homotopy type as Ext(S', R(e)). 
The claim now follows since one can extend any triangulation over R to a 
triangulation over R(e). □ 

Let C S'^ be the basic semi-algebraic set defined by 

S'^ = {x G R(e)'=+^ I ||x|| = 1, P/^(x) > 0,...,P^(x) > 0,P^+i(x) > 0}. 

Lemma 2.45. For < i < k, we have 

Proof. Note that Sb is bounded by Proposition 12.441 and S'^ is the 
projection from the origin of the set {1} x Sb C {1} x R(e)'^ onto the unit 
sphere S'^ in R(e)^''^^. Since Sb is bounded, the projection does not intersect 
the equator and consists of two disjoint copies (each homeomorphic to the 
set Sb) in the upper and lower hemispheres. □ 



CHAPTER 3 



Bounding the Betti Numbers 

1. Results 

We prove the following theorem. 

Theorem 3.1. [T7j Lei P = {Pi, . . .,Pm} C R[Xi, . . .,Xk], m < k. Let 
5 C be defined by 

Pl>0,...,Pm>0 

with deg(Pj) < 2. Then, forO<i<k-l, 

min{m.+l,k~i} . , .\ / , \ n / a 1 

As a consequence of Theorem 13. II we get a new bound on the sum of the 
Betti numbers, which we state for the sake of completeness. 

Corollary 3.2. Let V = {Pi, . . . ,Pm} C R[Xi, . . . ,Xk], m < k. Let 
S* C R^ be defined by 

Pl>0,...,Pm>0 

with deg(Pj) < 2. Then 

( min{m+l,k-i} . . .\ / , x\ 

Remark 3.3. The technique used in this chapter was proposed as a 
possible alternative method by Barvinok in [9], who did not pursue this 
further in that paper. Also, Benedetti, Loeser, and Risler |26j used a similar 
technique for proving upper bounds on the number of connected components 
of real algebraic sets in R'^ defined by polynomials of degrees bounded by 
d. However, these bounds (unlike the bounds we obtain) are exponential 
in k. Finally, there exists another possible method for bounding the Betti 
numbers of semi-algebraic sets defined by quadratic inequalities, using a 
spectral sequence argument due to Agrachev [1]. However, this method 
also produces a non-optimal bound of the form kP'<^^ (similar to Barvinok's 
bound) where the constant in the exponent is at least two. We omit the 
details of this argument referring the reader to |13j for an indication of the 
proof (where the case of computing, and as a result, bounding the Euler- 
Poincare characteristics of such sets is worked out in full details). 
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3. BOUNDING THE BETTI NUMBERS 



2. Proof Strategy 



Our strategy for proving Theorem 13.11 is as follows. Using certain in- 
finitesimal deformations we first reduce the problem to bounding the Betti 
numbers of another closed and bounded semi-algebraic set defined by a new 
family of quadratic polynomials. We then use inequalities obtained from the 
Mayer- Vietoris exact sequence to further reduce the problem of bounding 
the Betti numbers of this new semi-algebraic set to the problem of bound- 
ing the Betti numbers of the real projective algebraic sets defined by each 
£-tuple, £ < m, of the new polynomials. The new family of polynomials 
also has the property that the complex projective algebraic set defined by 
each ^-tuple, I < k, ol these polynomials is a non-singular complete inter- 
section. According to Theorem 12.411 we have precise information about the 
Betti numbers of these complex complete intersections. An application of 
the Smith inequality (see Theorem l2.22|) then allows us to obtain bounds on 
the Betti numbers of the real parts of these algebraic sets and, as a result, 
on the Betti numbers of the original semi-algebraic set. 

3. Constructing Non-singular Complete Intersections 

In Chapter 13.21 we introduced the notion of a projective complex alge- 
braic set which is a non-singular complete intersection (see Definition I2.39P . 
Next, we show the existence of such a set and how to obtain a non-singular 
complete intersection from a given algebraic set in complex projective space. 

Proposition 3.4. There exists a family H = {Hi, ... , Hm} C R[Xo, ... , X^] 
of positive definite quadratic forms such that Zer{Hj,FQ) is a non-singular 
complete intersection for every J C {1, . . . , m}. 

Proof. Recall that the set of positive definite quadratic forms is open 
in the set of quadratic forms over R. Moreover, any real closed field contains 
the real closure of Q. Thus, we can choose a family Ti. = {Hi, . . . , Hm} C 
R[Xo, . . . , Xk] of positive definite quadratic forms such that their coefficients 
are algebraically independent over Q. It follows by Bertini's Theorem (see 
[50] . Theorem 17.16) that Zer(Wj, J C {1, . . . is a non-singular 
complete intersection. □ 

The following proposition allows us to replace a family of real quadratic 
forms by another family obtained by infinitesimal perturbations of the orig- 
inal family and whose zero sets are non-singular complete intersections in 
complex projective space. 

Proposition 3.5. Let 



Q = {Qi,...,Qm}cR[Xo. 
be a set of quadratic forms and let 

n = {Hi,...,Hm} cR[Xo 




Xk] 
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be a family of positive definite quadratic forms such that Zer('H,P^) is a 
non-singular complete intersection for every J C {1, . . . , m}. 
Let 1 ^ 5 > be infinitesimals, and let 

Q = {Qi, ■ ■ ■ ,Qm} with 
Qi = {l-5)Qi + 5H,. 

Then for any J C {1, . . . , m}, 

ZeT{Qj,¥^C{5)) 
is a non-singular complete intersection. 

Proof. Consider 

Qt = {Qt,i, ■ ■ ■ ,Qt,m} with 
Qt,i = {l-t)Qi + tHi. 
Let J C {1, . . . , m}, and let Tj C C be defined by, 

Tj = {t G C I Zer(Qt^j, P(^) is a non-singular complete intersection }. 

Clearly, Tj contains 1. Moreover, since being a non-singular complete in- 
tersection is a stable condition, Tj must contain an open neighborhood of 
1 in C and so must T = Hj^-ji^ j^jT/. Finally, the set T is constructible, 
since it can be defined by a first order formula. Since a constructible sub- 
set of C is either finite or the complement of a finite set (see for instance, 
|19] . Corollary 1.25), T must contain an interval (0,to))*o > 0. Hence, its 
extension to C{5) contains 5. □ 

4. Proof of Theorem 13.11 

Before we prove Theorem 13.11 we need what follows next: 
Let P = {Pi, . . . , P„J C R[Xi, ...,Xk],m<k, with deg(Pi) < 2, 1 < i < 
m. Let S" C R*^ be the basic semi-algebraic set defined by 

5 = {x G R^^ I Pi(x) > 0, . . .,Pm{^) > 0}. 

Let 1 ^ e ^ 5 > he infinitesimals, and let 



Pm+l — I — £ y Xi . 



i=l 



Let Sh C R(e)'^ be the basic semi-algebraic set defined by 

5b = {xGR(e)'= I Pi(x) >0,...,P^(x) >0,P™+i(x) >0}. 

The homology groups of S and Sh are isomorphic by Proposition l2.44l More- 
over, the set Sh is bounded. 

Let C S*^ be the basic semi-algebraic set defined by 

S'^ = {x G R(e)'=+i I |x| = 1, Pi'^(x) > 0, . . . ,P^(x) > 0,P^+i(x) > 0}. 
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Then, for < i < k, we have 

6i(5fe,Z2) = hi{S\Z2). 

by Lemma 12.451 

We now fix a family of polynomials that will be useful in what fol- 
lows. By Proposition 13.41 we can choose a family Ti. = {Hi, . . . , Hm+i} C 
R[Xo, . . . ,Xk] of positive definite quadratic forms such that Zer('Hj, P^^^^) 
is a non-singular complete intersection for every J C {1, . . . ,m + 1}. 

Let Pi = (1 - 6)P!' + 6Hi, l<i<m + l. Let T (resp., f) be the basic 
semi-algebraic set defined by 

T = {xGR(e,5)'=+^ I ||x|| = 1, A(x) >0,...,P^(x) >0,,P^+i(x) >0} 
and 

T = {xGR(e,5)'=+' I ||x|| = 1, Pi(x) >0,...,P^(x) >0,P™+i(x) >0}, 

respectively. 
Also, let 

^ = {-Pl , ■ ■ ■ , Pm , Pm.+l } • 

Lemma 3.6. We have, 

(1) the homology groups of and T are isomorphic, 

(2) the homology groups of T and T are isomorphic, 

(3) for allJ C {!,..., m-M}, 

Zer(P7, P^^^ is a non-singular complete intersection, and 

(4) for allJ C {l,...,m + l}, 

hi (zer(pJ,Ext(S^R(e,5)),Z2) < 2h (zer(Pj, P^^^ 5^), Z2) . 

Proof. For the first part note that the sets Ext(S''', R(e, 5)) and f have 
the same homotopy type using Lemma 16.17 in |20| . 

The second part is clear since we have a retraction from T to T. 
The third part follows from Proposition 13.51 

For the last part, let tt : Ext(S'', R(e, 5)) "^Kie 5) double 
cover obtained by identifying antipodal points. Then the restriction of vr 
to Zer(:Pj,Ext(S*^,R(e,(5))) gives a double cover, 

TT : Zer(pJ,Ext(S^R(e,5))) ^ Zer(Pj, P^^,^^^). 

Now apply Proposition 12.24] □ 

Proposition 3.7. For < i < k — 1, we have 

min{m+l,fc— i} 

bi{T,Z2) <l + 2{k-m)+ 2^ 

j=0 



fm-\-l\f k \ 
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Proof. First note that by Lemma[3^ (j3|) Zer('Pj,P^^^ is a complete 
intersection for all J C {1, . . . , m + 1}. For < i < A; — 1, we have 



771 + 1 



bi{T,Z2) < 6i Ext(S^R(e,5))\ J Zer(P^,Ext(S^R(e,5))),Z2 

V i=l / 

/m+1 \ 

< l + ( U Zer(Pi,Ext(S^R(e,5))),Z2J , 



where the first inequality is a consequence of the fact that, T is an open 
subset of 



m+l 

Ext(S^R(e,5)) \ IJ Zer(P^,Ext(S^R(e,(5))) 

i=l 



and disconnected from its complement in Ext(S'', R(e, 5))\Urit^ '^^'^{Pi, Ext(S'', R{e 
and the last inequality follows from Theorem 12.231 (Alexander Duality) . 

It follows from Proposition l2 . 19 l (j2.2p . Lemma [3.6l (|^ and Theorem l2.22l (Smith 
inequality) that 



k—i 

bi{T,Z2) < 1 + 5^ 5^ 6fc_._,•(zer(pJ,Ext(S^R(e,5))),Z2) 
i=i \J\=j 

k—i 

i=i \J\=j 

m,in{m,+l,k—i} 

^1 + 2- E EK^''(^^''^cm>).^2). 
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Note that for j < m + 1 the number of possible j-ary intersections is equal 
to ("^^^) and using Theorem 12.411 we conclude 

min{m+l,k—i} , , 

6,(r,Z2) < 1 + 2- 

i=i V J / 

mm{m+l,fc— j} , , / j \ 

= l+2(fc + l)-2(m + l)+ 2\"DG-l) 

□ 

We are now in a position to prove Theorem 13.11 

Proof of Theorem 13.11 It follows from the Universal Coefficients The- 
orem (see |51j . Corollary 3.A6 (b)), that bi{S) < Z2). We have by 
Lemma 13.61 that the homology groups (with Z2 coefficients) of and T 
are isomorphic. Moreover 6j(5, Z2) = ^bi{S^ ,^2), for < i < A; — 1, by 
Proposition 12.441 and Lemma 12.451 Hence, the first inequality follows from 
Proposition 13.71 

The second inequality follows from an easy computation. □ 



CHAPTER 4 



Bounding the Number of Homotopy Types 

1. Result 

We prove the following theorem. 

Theorem 4.1. [16] Let R be a real closed field and let 

V = {Pi,...,Pm}cR[Yi,...,Ye,Xi,...,Xk], 

with degy(Pi) < 2,degx{Pi) < d,l < i < m. Let ir : R^+'= ^ R'^ be the 
projection on the last k co-ordinates. Then for any V-closed semi-algebraic 
set S C R^"*"*^, the number of stable homotopy types (see Definition \2.28\) 
amongst the fibers, S^, is bounded by {2"^ikd)'-^^"'''\ 

Remark 4.2. (1) The bound in Theorem 14.11 (unlike the one in 
Theorem I2.42p is polynomial in i for fixed m and k. The expo- 
nential dependence on m is unavoidable, as can be seen from a 
slight modification of Example 12.431 Consider the semi-algebaic 
set S C R^+^ defined by 

Yi{Yi - 1) = 0, l<i<m<e, 
0<X <Yi + 2-Y2 + ... + 2"^-! • Ym. 

Let vr : R^"*"^ — > R be the projection on the X-coordinate. Then, 
the sets S'x, x G {0, 1 . . . , 2"^"^}, have different number of connected 
components, and hence have distinct (stable) homotopy types. 
(2) The technique used to prove Theorem 12.421 in |22j does not di- 
rectly produce better bounds in the quadratic case, and hence we 
need a new approach to prove a substantially better bound in this 
case. For technical reasons, we only obtain a bound on the number 
of stable homotopy types, rather than homotopy types. But note 
that the notions of homeomorphism type, homotopy type and sta- 
ble homotopy type are each strictly weaker than the previous one, 
since two semi-algebraic sets might be stable homotopy equivalent, 
without being homotopy equivalent (see [73j . p. 462), and also ho- 
motopy equivalent without being homeomorphic. However, two 
closed and bounded semi-algebraic sets which are stable homotopy 
equivalent have isomorphic homology groups. 
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4. BOUNDING THE NUMBER OF HOMOTOPY TYPES 



2. Proof Strategy 

The strategy underlying our proof of Theorem 14.11 is as follows. We 
first consider the special case of a semi-algebraic subset, A C S^, defined 
by a disjunction of m homogeneous quadratic inequalities restricted to the 
unit sphere in R^"*"^. We then show that there exists a closed and bounded 
semi-algebraic set C (see (|4.14p below for the precise definition of the semi- 
algebraic set C"), consisting of certain sphere bundles, glued along certain 
sub-sphere bundles, which is homotopy equivalent to A. The number of these 
sphere bundles, as well descriptions of their bases, are bounded polynomially 
in i (for fixed m). 

In the presence of parameters Xi, . . . , X^, the set A, as well as C , will 
depend on the values of the parameters. However, using some basic homo- 
topy properties of bundles, we show that the homotopy type of the set C 
stays invariant under continuous deformation of the bases of the different 
sphere bundles which constitute C. These bases also depend on the param- 
eters, Xi, . . . , Xk, but the degrees in Xi, . . . , X^ of the polynomials defining 
them are bounded by 0{id). Now, using techniques similar to those used in 
[22], we are able to control the number of isotopy types of the bases which 
occur as the parameters vary over R'^. The bound on the number of isotopy 
types, also gives a bound on the number of possible homotopy types of the 
set C", and hence of A, for different values of the parameter. 

In order to prove the results for semi-algebraic sets defined by more gen- 
eral formulas than disjunctions of weak inequalities, we first use Spanier- 
Whitehead duality to obtain a bound in the case of conjunctions, and then 
use the construction of homotopy colimits to prove the theorem for gen- 
eral 7^-closed sets. Because of the use of Spanier- Whitehead duality we get 
bounds on the number of stable homotopy types, rather than homotopy 
types. 



3. Topology of Sets Defined by Quadratic Constraints 

One of the main ideas behind our proof of Theorem 14. II is to parametrize 
a construction introduced by Agrachev in [l] while studying the topology of 
sets defined by (purely) quadratic inequalities (that is without the parame- 
ters Xi, . . . ,Xk in our notation). However, we avoid construction of Leray 
spectral sequences as was done in [l]. For the rest of this section, we fix a 
set of polynomials 

Q = {Qi,...,Qm}cR[Yo,...,Y,,Xi,...,Xk] 

which are homogeneous of degree 2 in Yq, . . . , 1^, and of degree at most d in 
Xi, . . . , Xk. 

We will denote by 

Q=(Qi,...,Qm):R'+'xR^^R-, 
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the map defined by the polynomials Qi, ■ ■ ■ ,Qm, and generally, for / C 
{1, . . . , m}, we denote by Qi : R^+^ x R'' ^ R^, the map whose co-ordinates 
are given hy Qi, i £ I. When / = [m], we will often drop the subscript I 
from our notation. 

For any subset I C [m], let Aj G x be the semi-algebraic set 
defined by 

(4.1) Ai = |J{(y,x) I |y| = 1 A Qi(y,x) < 0}, 

iei 

and let 

(4.2) 0/ = {w G R™ I |u;| = l,uji = 0,i ^ 1 ,0Ji < 0,i e I}. 

For to £ Qj we denote by ljQ G R[1o; ■ ■ ■ ,Ye, Xi, . . . , X^] the polynomial 
defined by 

m 

(4.3) ojQ = ^uJiQi. 

i=0 

For (u;,x) S F/ = 0/ x R''^ we will denote by ujQ{-,x) the quadratic 
form in Yq, . . . ,Yi obtained from loQ by specializing = Xj, 1 < i < /c. 
Let Bi C nj X S'^ X R^ be the semi-algebraic set defined by 

(4.4) Bi = {(a;, y, x) I (J G ^i, y G S^ x G R^ u;Q{y, x) > 0}. 

We denote by <f>i : Bj ^ Fj and (f)2 ■ Bj ^ x r!' the two projection 
maps (see diagram below). 




The following key proposition was proved by Agrachev [T] in the unparametrized 
situation, but as we see below it works in the parametrized case as well. 

Proposition 4.3. The map 02 gives a homotopy equivalence between 
Bj and (l)2{Bi) = Aj. 

Proof. In order to simplify notation we prove it in the case I = [m], 
and the case for any other / would follow immediately. We first prove that 
4>2{B) = A. If (y,x) G A, then there exists some i,l < i < m, such that 
Qi(y, x) < 0. Then for uj = (— . . . , —6m,i) (where 5ij = 1 if i = j, and 
otherwise), we see that (ti;,y,x) G B. Conversely, if (y,x) G (j)2{B), then 
there exists uj = {toi, . . . ,uim) G ^ such that, X^^i i^iQi(y, x) > 0. Since 
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"^i < 0, 1 < i < m, and not all LOi = 0, this implies that (5j(y,x) < for 
some i,l < i < m. This shows that (y, x) G A. 
For (y, x) G 02 the fiber 

(^^""^(y, x) = {{uj, y, x) I a; G such that ojQ{y, x) > 0} 

is a non-empty subset of O defined by a single linear inequality. Thus each 
non-empty fiber is an intersection of a convex cone with S™""^, and hence 
contractible. 

The proposition now follows from the well-known Vietoris-Smale theo- 
rem (see Theorem I2.27p . □ 

We will use the following notation. 

Notation 4.4. For any quadratic form Q G R[yoi • • • , we will denote 
by index((5) the number of negative eigenvalues of the symmetric matrix of 
the corresponding bilinear form, that is of the matrix Mq such that, Q{'y) = 
(Mgy, y) for all y G R^"*"^ (here (•, •) denotes the usual inner product). We 
will also denote by \i{Q), < i < i, the eigenvalues of Q in non-decreasing 
order, i.e., 

Ao(Q)<Ai(Q)<---<AKQ). 

For / C [m], let 

(4.6) Fij = {(a;,x) e fti x R*^ | index(a;Q(-, x)) < j}. 

It is clear that each Fj j is a closed semi-algebraic subset of Fj and that 
they induce a filtration of the space Fj given by 

Fifl C C ••• C Fi^i+i = Fi. 

Lemma 4.5. The fiber of the map cpj^i over a point {uj, x) G Fjj \ Fjj-i 
has the homotopy type of a sphere of dimension i — j. 

Proof. As before, we prove the lemma only for / = [m]. The proof for 
a general / is identical. First notice that for (w,x) G Fj \ Fj^i, the first j 
eigenvalues of loQ{-,x) 

Ao(a;Q(-,x)),... ,Aj_i(a;Q(-,x)) < 0. 

Moreover, letting Wo{ujQ{-,x)), . . . ,Wi{ioQ{-,x)) be the co-ordinates with 
respect to an orthonormal basis eo{ujQ{-,x)), . . . , ei{LoQ{-,x)), consisting of 
eigenvectors of ujQ{-,x), we have that (f)i^{uj,x) is the subset of = {uj} x 
X {x} defined by 

e 

^Ai(a;Q(-,x))W,,(a;Q(-,x))2 >0, 

i=0 

e 

J^Ty,(c^Q(-,x))2 = l. 

1=0 
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Since, Xi{LoQ{-,x)) < 0, < i < j, it follows that for (u;,x) G Fj \ Fj-i, 
the fiber cf)^^ {oj , :x.) is homotopy equivalent to the — j)-dimensional sphere 
defined by setting 

WoiioQi;^)) = ■■■ = VF,_i(a;Q(-,x)) = 
on the sphere defined by J2i=o ^ii^Qi'i'^))'^ = 1- 1^ 

For each (ti;,x) G Fjj \ let L^(u;,x) C R^'^^ denote the sum of 

the non-negative eigenspaces of ujQ{-,:x.) (i.e., L^(u;,x) is the largest lin- 
ear subspace of R^"*"^ on which ujQ{-,x) is positive semi-definite). Since 
index(Li;(5(") x)) = j stays invariant as {u>, x) varies over Fjj\Fjj^i, L^i^^, x) 
varies continuously with (u;,x). 

We will denote by Cj the semi-algebraic set defined by 

i+i 

(4.7) Ci = |J{Ky>x) I i^^^) e Fij\Fi^,^i,yG L+(c^,x),|y| = 1}. 
j=0 

The following proposition relates the homotopy type of Bj to that of C/. 

Proposition 4.6. The semi- algebraic set Cj defined above is homotopy 
equivalent to Bj (see lli4-4\ ) f'^^ definition of Bj). 

Proof. We give a deformation retraction of Bj to Cj constructed as 
follows. For each (cj, x) £ Fjj \ we can retract the fiber (j)'^'^{uj, x) to 

the zero-dimensional sphere, L^{uj,x) H by the following retraction. Let 

Wo{ujQl{;x)),...,Wi{LoQl{;x)) 

be the co-ordinates with respect to an orthonormal basis eo{ujQ{-,x)), . . . , e^(u;(5(-, x)), 
consisting of eigenvectors of ujQi{-,x) corresponding to non-decreasing order 
of the eigenvalues of ujQ{-,x). Then, (f)^^{oj,x) is the subset of defined 
by 

e 

^Xi{ujQl{;x))W^{L0Ql{;x)f > 0, 

£ 

J2Wiiu;Qi{;x)f = 1. 

i=0 

and L~^{u>,x) is defined by Wq{ujQi{-,x)) = ••• = Wi-\{ujQi{-,x)) = 0. 
We retract (f)'^^{uj^x) to the zero-dimensional sphere, L^{uj,x) D by the 
retraction sending, 

{wo, ...,wi)e (p^^{uj,x), 

at time t to 

((1 - t)wo, . . . , (1 - t)wi-i,t'wi), 
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where < t < 1, and 



t' 



i-(i-i)'E: 



i-i 

i=0 



1/2 



Notice that even though the local co-ordinates {Wq, . . . , Wi) in R^"*"^ with 
respect to the orthonormal basis (cq, • • • , e^) may not be uniquely defined at 
the point {io,x) (for instance, if the quadratic form ujQi{-,x) has multiple 
eigenvalues), the retraction is still well-defined since it only depends on the 
decomposition of into orthogonal complements span(eo, • • • , e£„i) and 
span(e£). We can thus retract simultaneously all fibers over Fji \ Fj^i-i 
continuously, to obtain a semi-algebraic set Bj £ C Bj, which is moreover 
homotopy equivalent to Bj. 

This retraction is schematically shown in Figure [U where Fj^i is the 
closed segment, and -F/,^-! are its end points. 



Figure 1 . Schematic picture of the retraction of Bj to Bj , 



Now starting from i?//, retract all fibers over Fj^£-i \ -F//_2 to the 
corresponding one dimensional spheres, by the retraction sending 

at time t to 

((1 - t)wo, . . . , (1 - t)wi_2,t'we-i,t'we), 
where < i < 1, and 



E 



i=i-i 
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to obtain which is homotopy equivalent to i?//. Continuing this 

process we finally obtain Bj^ = Ci, which is clearly homotopy equivalent 
to Bj by construction. □ 

Notice that the semi-algebraic set (f)^^{Fij \ Fjj^i) n C/ is a S^"-'- 
bundle over Fij \ Fjj-i under the map (pi, and Cj is a union of these 
sphere bundles. We have good control over the bases, Fjj \ of these 

bundles, that is we have good bounds on the number as well as the degrees of 
polynomials used to define them. However, these bundles could be possibly 
glued to each other in complicated ways, and it is not immediate how to 
control this glueing data, since different types of glueing could give rise to 
different homotopy types of the underlying space. In order to get around this 
difficulty, we consider certain closed subsets, Fj ■ of F/, where each F[ ■ is an 

infinitesimal deformation of Fij\Fij-i, and form the base of a S^"-' -bundle. 
Moreover, these new sphere bundles are glued to each other along sphere 
bundles over Fj j n Fj -^^^ and their union, C|, is homotopy equivalent to 
C/. Finally, the polynomials defining the sets F[ ■ are in general position in 
a very strong sense, and this property is used later to bound the number of 
isotopy classes of the sets Fl ■ in the parametrized situation. 

We now make precise the argument outlined above. Let A/ be the poly- 
nomial in R[Zi, . . . , Zm, Xi, . . . , X^., T] defined by 

A/ = det(Mz,.o + rid£+i), 

where Zj ■ Q = ^-^j ZiQi, and each Hjj G R[Zi, ... , Zm, Xi,. . . , Xk]. 

Notice, that Hj j is obtained from Hj = H^^j j by setting the variable 
Zi to in the polynomial Hj for each i ^ L 

Note also that for (z,x) G R™ x R*^, the polynomial A7(z,x, T) being 
the characteristic polynomial of a real symmetric matrix has all its roots 
real. It then follows from Descartes' rule of signs (see for instance j20j ). 
that for each (z,x) G R"^ x R*^, where Zj = for all i ^ I, index(z(5(-, x)) is 
determined by the sign vector 

(sign(i?/_^(z, x)), . . . , sign(i?/,o(z, x))). 

Hence, denoting by 

(4.8) TCi = {Hi ft, . . . C R[Zi, . . . ,Zm,Xi, . . . ,Xk], 
we have 

Lemma 4.7. For each j,0 < j < i + l, Ffj is the intersection of Fj with 
a Til-closed semi- algebraic set Dj j C R'""'''^. 

Notation 4.8. Let Djj be defined by the formula 

(4.9) Dij= U 
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for some T,jj C Sign('H/). Note that, Sign{7ii) C Sign('H) and T,jj C 
for all I C [m]. 

Now, let = {6i, . . . ,5o) and e = (e^+i, . . . , eo) be infinitesimals such 
that 



Given a G Sign(?^/), and < j < £ + 1, we denote by 7^(^Tp C R'™+^ 
the set defined by the formula aj obtained by taking the conjunction of 

—Ej — 5i < Hj^i < £j + 5i for each Hj^i G TLj such that a{Hj^i) = 0, 
Hj^i > —Ej — Si, for each H[ i G Tij such that a{Hj^i) = 1, 
Hi^i < Ej + 6i, for each Hj i G 7ii such that a{Hj^i) = —1. 

Similarly, we denote by C R'™+'= the set defined by the formula 

(7° obtained by taking the conjunction of 

—Ej — 5i < Hji < Ej + 5i for each Hij G 7ii such that a{Hj^i) = 0, 
Hj i > —Ej — 5i, for each Hj^i G such that a{Hf^i) = 1, 
< + for each Hf i G such that a{Hi^i) = —1. 

For each j, < j < ^ + 1, let 



< (^0 < • • • < < eo < • • • < ££+1 < 1 



and let 



(4.10) 



R' = R(e, 5) 



Dh = U ^(^|)' 




Lemma 4.9. For eachj,0 < j < i + l, D'j ■ is aTC'j-closed semi-algebraic 
set, where 



e e+i 



(4.12) 




i=Ojr=0 



Proof. Follows from the definition of the sets D'j 
Lemma 4.10. For 0<j + l<i<£ + l, 



□ 



^;,.nz?;,, = 0. 
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Proof. In order to keep notation simple we prove the proposition only 
for / = [m]. The proof for a general / is identical. The inclusions, 

Dj^i C Dj c A-i c A, 
follow directly from the definitions of the sets 
and the fact that, 

£j-l < Ej < Ei-i < £i. 

It follows immediately that, 

D[ = Di\DU 

is disjoint from Z)^, and hence from D'-. □ 

We now associate to each F'j ■ a. {I — j)-dimensional sphere bundle as 
follows. For each (a;,x) G F^ - = Fjj \ F'j -^-^, let L+(a;,x) C R^+^ denote 
the sum of the non-negative eigenspaces of ujQ{-,'x) (i.e., L'j{uj,x.) is the 
largest linear subspace of R^"*"^ on which w(5(-,x) is positive semi-definite). 
Since mdex{ujQ{-,:x.)) = j stays invariant as (w,x) varies over F'/^, L^{lj,:x.) 
varies continuously with (t(j,x). 

Let, 

Ao(ti-', x) < • • • < Aj_i(w, x) < < Aj(ci;, x) < • • • < Xi{uj, x), 

be the eigenvalues of ujQ{-,x.) for (i^,x) G ^i'j- There is a continuous ex- 
tension of the map sending (u;,x) h->- Lt(a;,x) to (w,x) S Fj j. 

To see this observe that for (uj, x) G F'/ j the block of the first j (neg- 
ative) eigenvalues, Ao(w,x) < ••• < Aj_i(a;,x), and hence the sum of the 
eigenspaces corresponding to them can be extended continuously to any in- 
finitesimal neighborhood of Fjj, and in particular to F'j j. Now Lj'{uj,x.) is 
the orthogonal complement of the sum of the eigenspaces corresponding to 
the block of negative eigenvalues, Ao(u;,x) < • • • < Aj_i(lj,x). 

We win denote by C'j j C F'j j x R'^+^ the semi-algebraic set defined by 

(4.13) C'jj = {{to, y, x) I {oj, x) G F'j^^,y G L+{oj, x), |y| = 1}. 

Note that the projection ttjj : C'j j ^ Fj makes C'jj the total space 
of a (^ — j)-dimensional sphere bundle over F'^ j. 
Now observe that 

c'jj_,nc'j^^=7rY]{FijnF'j^^_,), 

and 

is also a (£ — j) dimensional sphere bundle over Fj j n F[ 
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Let 

(4.14) = U C'l^^. 

j=0 

We have that 

Proposition 4.11. C| is homotopy equivalent to Ext(C/,R'), where Cj 
and R' are defined in \4- 'A ) ^.i^d J^.iQp respectively. 

Proof. Let e = (e^+i, . . . , eo) and let 

' R{e,6i,...,5i), 0<i<i, 
Ri = \ R{ei+i, Ei-t^i), i + l<i<2£ + 2, 
R,i = 2i + 3. 

First observe that Cj = hmg^^j Cj where Cj is the semi-algebraic set defined 
in (j4.7p above. 
Now let, 

= C'l, 
Cifl = limCj, 

Cj^i = limC/.j_i, 1 < i < 

Si 

C/,i = \im Ci^i-i,l + 2<i<2i + 3. 

Notice that each Cj^i is a closed and bounded semi-algebraic set. Also, 
for i > 0, let C/,i_i,i C R7^+^+*^ be the semi-algebraic set obtained by 
replacing 5i (resp., ej) in the definition of C/^j_i by the variable t. Then, 
there exists to > 0, such that for all < ti < t2 < to, Cj^i-i^n C C/^.j_i^t2- 

It follows (see Lemma 16.17 in |20| ) that for each i, < i < 2£ + 3, 
Ext(C/^i,R,) is homotopy equivalent to C/.i_i. □ 

4. Partitioning the Parameter Space 

The goal of this section is to prove the following proposition (Proposition 
I4.12P . The techniques used in the proof are similar to those used in |22j for 
proving a similar result. We go through the proof in detail in order to extract 
the right bound in terms of the parameters d,k,i and m. 

Proposition 4.12. There exists a finite set of points T C R''' with 

#r < (2™£A;d)^('"*^) 

such that for any x G R*^, there exists z G T, with the following property. 

There is a semi- algebraic path, 7 : [0, 1] — > R''^ and a continuous semi- 
algebraic map, (j) : ^} X [0,1] ^ il. (see j^.S^ and ( f^. for the definition of 
and R' 
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with 7(0) = X, 7(1) = z, and for each I C [m], 

is a homeomorphism for each < t < 1. 

Before proving Proposition 14.121 we need a few preliminary results. Let 

(4.15) n" = n'u{Zu..., Zm., zl + --- + zl- 1}, 

where 7i' = is defined in ()4.12p above. 

Note that for each j, < j < Fj j is a. W'-closed semi- algebraic set. 
Moreover, let i/j : — > Ji'^ be the projection onto the last k co-ordinates. 

Notation 4.13. We fix a finite set of points T C R'' such that for every 
X € R'^ there exists z G T such that for every "^''-semi-algebraic set V, the 
set ^p^^{x.) n y is homeomorphic to ijj^^{z) H V. 

The existence of a finite set T with this property follows from Hardt's 
triviality theorem (Theorem 12. 15 p and the Tarski-Seidenberg transfer prin- 
ciple, as well as the fact that the number of "^''-semi-algebraic sets is finite. 

Now, we note some extra properties of the family 7i" . The notations 
SigUp and TZ{a) were introduced in Chapter 11.11 

Lemma 4.14. If a £ Signp(W"), then p < k + m and n{a) C R'™+*^ is 
a non-singular (m + k — p)- dimensional manifold such that at every point 
(z,x) G 'R-{a), the (p x (m -|- k))-Jacobi matrix, 

f dP dP \ 

\dZi dYj J p^-j^i,^ cr(P)=0, l<i<m., l<j<fc 

has maximal rank p. 

Proof. Let Ext(S''"~^,R') be the unit sphere in i?"". Suppose without 
loss of generality that 

m 

{P G n"\ a{P) = 0} = {H,, - e,, - 6^,, . . . , H,^_, - e,^_, - S,^_, , J] - 1} 

1=1 

since the equation Zi = eliminates the variable Zi from the polynomials. 
It follows that it suffices to show that the algebraic set 
p-i 

(4.16) V = n{(z,x) G Ext(S™-\R') X R"= | i/,,(z,x) = Sj^ + 6,^} 

r=l 

is a smooth ((m — 1) + k — {p— l))-dimensional manifold such that at every 
point on it the {p x {m + A:))-Jacobi matrix, 

dP dP' 



dZi dYj y pg^//^ cr(P)=0, l<i<m, l<i<fc 



has maximal rank p. 
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Let p < m+k. Consider the semi-algebraic map : S™~^ xR'^ ^ 

RP-i defined by 

(z,x) ^ (Fi,(z,x), . . . ,iJip_^(z,x)). 

By the semi-algebraic version of Sard's theorem (see |27j ). the set of critical 
values of -Pji,...,ip_i is a semi-algebraic subset C of R^~^ of dimension strictly 
less than p — 1. Since 5 and e are infinitesimals, it follows that 

{sj, +6i„..., ej^_, + 6,^_,) i Ext(C, i?')- 

Hence, the algebraic set V defined in (I4.16P has the desired properties, and 
the same is true for the basic semi-algebraic set TZ{a). 

We now prove that p < m + k. Suppose that p > m + k. As we have 
just proved, 

{i7i,(z,x) =ej, +5i,,...,if,„^^,_,(z,x) = ei„+fe_i + J 

is a finite set of points. But the polynomial Hip_-^ — £jp_i — cannot 
vanish on each of these points as 6 and e are infinitesimals. □ 

Lemma 4.15. For every x G R^, and a £ Sigiip{TC'^) , where 

n'; = {p{Zi,...,Zm,^)\pen"}, 

the following holds. 

(1) < p < m, and 

(2) TZ{a)n'ilj~^{x.) is a non-singular (m — p)- dimensional manifold such 
that at every point (z,x) G '^(f) H '0~^(x), the {p x m)-Jacobi 
matrix, 

'dP\ 

has maximal rank p. 

Proof. Note that P^ = P{Zi, Zm, x) G R'[Zi, Z^] for each P £ 
n" andx G R^ The proof is now identical to the proof of Lemma I4.14[ □ 

Lemma 4.16. For any bounded 7i" -semi-algebraic set V defined by 
V= U ^(^)' 



CTeSvCSign(W") 



the partitions 



(TeSign(H") 



/m+k 



are compatible Whitney stratifications of R and V respectively. 

Proof. Follows directly from the definition of Whitney stratification 
(see [46l [34|), and Lemma HH □ 
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Fix some sign condition a S Sign('H"). Recall that (z,x) G 'R-{a) is a 
critical point of the map ip-jii^a) if the Jacobi matrix, 

(— ) 

at (z,x) is not of the maximal possible rank. The projection ijj{z,x.) of a 
critical point is a critical value of ipTi{u) ■ 

Let Ci C R'''"+'^ be the set of critical points of ip-jzia) over all sign con- 
ditions 

a G U Signp(W"), 

(i.e., over all a G Signp('H") with dim(7?.((T)) > k). For a bounded W-semi- 
algebraic set V, let Ci{V) C y be the set of critical points of ipnia) 0"^^^ 
sign conditions 

a G IJ Signp(7t:") n 

(i.e., over all a G with d.\m{JZ{a)) > k). 

Let C2 C R'"'+^ be the union of TZ{a) over ah 

0- G IJ Signp(W") 

(i.e., over all a G Signp('H") with dim(7^((T)) < k). For a bounded "^''-semi- 
algebraic set V , let C2{V) C F be the union of 7^(cr) over all 

G y Signp(?t:") n T.V 

p>m 

(i.e., over all a G with dim{Tl{a)) < k). 

Denote C = C7i U C2, and C7(y) = C7i(y) U C2{V). 

Lemma 4.17. For each hounded TC" -semi-algebraic V, the set C{V) is 
closed and bounded. 

Proof. The set C{V) is bounded since V is bounded. The union C2(y) 
of strata of dimensions less than k is closed since V is closed. 

Let (Ti G Signp^(W')nSv, <T2 G Sign^^ {H")r]T,v, where pi <m,pi< p2, 
and if <7i(P) = 0, then (T2{P) = for any P G H" . It follows that stratum 
7l{a2) lies in the closure of the stratum Tl{ai). Let ^7 be the finite family of 
{pi X pi)-minors such that Zer(j7',R') r\TZ{ai) is the set of all critical points 
of 7r7^((7i)- Then Zer(J',R') n 7^(o"2) is either contained in C2{V) (when 
dim(7^((72)) < k), or is contained in the set of all critical points of T^-jii^rj^) 
(when dim(7^((72)) > k). It follows that the closure of Zer(^7,R') n TZ{ai) 
lies in the union of the following sets: 

(1) Zer(J,R')n7^(al), 

(2) sets of critical points of some strata of dimensions less than m + 
k-pi, 
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(3) some strata of dimension less than k. 
Using induction on descending dimensions in case ([2]), we conclude that 
the closure of Zev{J',R') n 7^((Ti) is contained in C{V). Hence, C{V) is 
closed. □ 

Definition 4.18. We denote by Gi = il){Ci),i = 1, 2, and G = Gi U G2. 
Similarly, for each bounded "^''-semi-algebraic set V ^ we denote by GiiV) = 
^{Ci{V)), i = l,2, and G(y) = Gi{V) U G2iV). 

Lemma 4.19. We have TnG = 0. In particular, TnG{V) = for every 
bounded 7i" -semi-alqebraic set V . 



Proof. By Lemma H?T^ for all x G T, and a G Signp(7i^), 

(1) < p < m, and 

(2) 7^(cr) n V'~^(x) is a non-singular (m — p)-dimensional manifold such 
that at every point (z,x) G '^(o") H the {p x ?Ti)-Jacobi 
matrix, 

(—) 

has the maximal rank p. 
If a point xG TnGi = Tn ■(/'(Ci), then there exists z G R"" such that 
(z,x) is a critical point of '0-r,{o-) for some a G Up<m ^i^'i this is 

impossible by 

Similarly, x G TnG2 = Triip{C2), implies that there exists z G R'™ such 
that (z,x) G 7^(0") for some a G Up>m Signp('H"), and this is impossible by 
©• □ 

Let D be a connected component of R''^ \G, and for a bounded "^''-semi- 
algebraic set V, let D{V) be a connected component of ip^V) \ G{V). 

Lemma 4.20. For every bounded TC" -semi-algebraic set V, all fibers 
ip~^{x.) n y, X G -D are homeomorphic. 

Proof. Lemma [4.151 and Lemma [4.161 imply that V = i/j'^ {'ipiV) \ 
G{V)) n F is a Whitney stratified set having strata of dimensions at least 
k. Moreover, iply is a proper stratified submersion. By Thom's first isotopy 
lemma (in the semi-algebraic version, over real closed fields [34]) the map 
tpl^ is a locally trivial fibration. In particular, all fibers 'tp~^{'K) n y, x G 
D{V) are homeomorphic for every connected component D{V). The lemma 
follows, since the inclusion G{V) C G implies that either D C D{V) for 
some connected component D{V), or D n ip{V) = 0. □ 

Lemma 4.21. For each x G T, there exists a connected component D of 
R'^ \ G, such that ■(/)~^(x) DV is homeomorphic to D V for every 

bounded TC" -semi-algebraic set V and for every xi G D. 

Proof. Let y be a bounded "^''-semi-algebraic set and x G T. By 
Lemma [4.19t x belongs to some connected component D of R''^\G. Lemma [4.20l 
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implies that tp ^(x) n y is homeomorphic to iJj ^(xi) n V for every xi S 
D. □ 



We now are able to prove Proposition 14.12] 

Proof of Proposition 14.121 Recall that G = GiU G2, where Gi is 
the union of sets of critical values of ip-ji{u) over all strata 7^((t) of dimensions 
at least /c, and G2 is the union of projections of all strata of dimensions less 
than k. 

By Lemma r4.21l it suffices to bound the number of connected components 
of the set R''^ \ G. Denote by Ei the family of closed sets of critical points 
of ipz{a)j over all sign conditions a such that strata have dimensions 

at least k (the notation Z{a) was introduced in Chapter II. ip . Let £2 be the 
family of closed sets 2{a), over all sign conditions a such that strata TZ{a) 
have dimensions equal to A: — 1. Let £ = £1 U £2- Denote by E the image 
under the projection ip of the union of all sets in the family £. 

Because of the transversality condition, every stratum of the stratifica- 
tion of V, having the dimension less than m + k, lies in the closure of a 
stratum having the next higher dimension. In particular, this is true for 
strata of dimensions less than k — 1. It follows that G C E, and thus every 
connected component of the complement R* \ £' is contained in a connected 
component of R''^ \ G. Since dim(ii^) < k, every connected component of 
R"^ \ G contains a connected component of R''^ \ E. Therefore, it is suffi- 
cient to estimate from above the Betti number bo(R"^ \ E) which is equal to 
hk^i{E) by the Alexander's duality. 

The total number of sets Z{a), such that a G Sign('H") and dini{Z{a)) > 
k — 1, is 0{i'^^"^^^^) because each 2{a) is defined by a conjunction of at most 
m + 1 of possible 0{i'^ + m) polynomial equations. 

Thus, the cardinality il^£, as well as the number of images under the pro- 
jection TT of sets in £ is 0(£^(™'+^)). According to (12. 2p in Proposition 12.191 
hk-i{E) does not exceed the sum of certain Betti numbers of sets of the 
type 

n ^(^*)' 

l<i<p 

where every Ui € £ and 1 < p < k. More precisely, we have 

hk-i{E) < Y E bfc-M n ^(^^)) • 

i<p<fc {Ui,...,Up}c e \i<i<p ) 

Obviously, there are 0(^^(™+^)'') sets of the kind 

Using inequality ()2.3p in Proposition 12.191 we have that for each <^ as 
above, the Betti number bfc_p(<I>) does not exceed the sum of certain Betti 
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numbers of unions of the kind, 

^= U <U,^)=^i U 

i<j<q \l<j<q 

with 1 < q < p. More precisely, 

1<9<P l<n<-<i9<p \ \l<j<q 

It is clear that there are at most 2^ < 2^ sets of the kind ^. 

If a set U € £i, then it is defined by m polynomials of degrees at most 
0{£d). If a set [/ G £2, then it is defined by 0(2™) polynomials of degrees 
O^mid), since the critical points on strata of dimensions at least k are defined 
by 0(2™) determinantal equations, the corresponding matrices have orders 
0(111), and the entries of these matrices are polynomials of degrees at most 
olid). 

It follows that the closed and bounded set 

U 

i<j<q 

is defined by 0(A;2'")) polynomials of degrees 0{id). 

By Proposition [2IMI ^k-p+g-ii'^') < (2™fc^d)'^("''=) for all 1 < p < k, 
l<q<P- Then bfc_p(^>) < {I'^kid)^'^'"^^^ for every l<p<k. Since there 
are we get the claimed bound 

bfc_i(^) < {2'^kidf^'^''l 
The rest of the proof follows from Proposition 14.211 □ 



5. Proof of the Result 

5.1. The Homogeneous Case. We first consider the case where all 
the polynomials in Q are homogeneous in variables Yq, . . . ,Y£ and we bound 
the number of homotopy types among the fibers 5x, defined by the Q-closed 
semi-algebraic subsets 5 of x R'^. We first the prove the following theorems 
for the special cases of unions and intersections. 

Theorem 4.22. Let R be a real closed field and let 

Q = {Qi,...,Qm}cR[Yo,...,Yi,Xi,...,Xk], 

where each Qi is homogeneous of degree 2 in the variables Yq, . . . ,Y£, and of 
degree at most d in Xi , . . . , X^ ■ 

For i G [m], let C x R^' he semi- algebraic sets defined by 

^i = {(y,x) I |y| = 1 A Q,(y,x) <0)}, 

Let vr : X R'^ — > R'^ be the projection on the last k co-ordinates. 
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m 

Then, the number ofhomotopy types amongst the fibers |^ Aj^x is bounded 

1=1 

by 

With the same assumptions as in Theorem 14.221 we have 

Theorem 4.23. The number of stable homotopy types amongst the fibers 

m 

Aj,x is bounded by 

1=1 

(2"^ttd)°('"^). 

Before proving Theorems 14.221 and 14.231 we first prove two prehminary 
lemmas. 

Lemma 4.24. There exists a finite set T C R*^, with 

#r < (2'"£fc(i)°('"'=), 

such that for every x G R'^' there exists z G T, a semi- algebraic set Dx,z C 
R""''"^, and semi- algebraic maps /x,/z) o.s shown in the diagram below, such 
that /x, fz <ire both homotopy equivalences. 



(4.17) 




Ext( IJ Ai,x,R') Ext( U 

i£[m] i£[m] 

Moreover, for each I C [m], there exists a subset Dj^^^ C -Dx.z; such 
that the restrictions, /7,xi//,zi o/ /x, /z to Di,yL,z give rise to the following 
diagram in which all maps are again homotopy equivalences. 



-D/,x,z 

/l,z 

Ext(|J Ai,x, R') Ext(U Ai,z, R') 

For each / C J C [m], I?/,x,z C D j^x,z o-'ri'd the maps //^x, //,z a'^e restrictions 

0//j,x,/j,z- 

Proof of Lemma I4.24[ By Proposition l4.12l there exists T C R'^ with 

#r< {2'^ikdf^'^^\ 

such that for every x G R'^, there exists z G T, with the fohowing property. 
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There is a semi-algebraic path, 7 : [0, 1] R' and a continuous semi- 
algebraic map, (/) : X [0, 1] — > with 7(0) = x, 7(1) = z, and for each 
/ C [m], 

is a homeomorphism for each < t < 1 (see (|4.2p . (j4.10p and (j4.1ip for the 
definition of 0, R' and Fj j). 

Now, observe that Cjj^ (resp. Cjj^) is a sphere bundle over Fjj^ 
(resp. Fjj^). Moreover 

C/,i,x = {(^^^y) I ^ e -^/j-.x^y e L+(u;,x), |y| = 1}, 

and, for uj G Fjj^x ^ -^/j-i,x' '^^ ^^-^^ Lj'{ll),:x.) C L^_^(liJ,x). 

We now prove that the map (p induces a homeomorphism (j) : C'^ — > 
C^, which for each / C [m] and < j < £ restricts to a homeomorphism 

^id ■ ^i,j,z- 

First recall that by a standard result in the theory of bundles (see for 
instance, |42| . p. 313, Lemma 5), the isomorphism class of the sphere bundle 
Cj-^ Ft ■, is determined by the homotopy class of the map, 

Flj,. ^ Gr{£ + l-j,£ + l) 
UJ ^ L+(w,x), 

where Gr{m,n) denotes the Grassmannian variety of m dimensional sub- 
spaces of R'". 

The map (p induces for each j, < j < £, a homotopy between the maps 

fo-F'jj^^ ^ Gr(^ + l-j,^ + l) 
Lu 1— > L^(w,x) 

and 

^ Gr{£+l-j,£+l) 
LO ^ Lj'{lo,z) 

(after indentifying the sets Fj ■ ^ and ■ ^ since they are homeomorphic) 
which respects the inclusions L^(a;, x) C L^_]^(a;,x), andL^(u;,z) C L^_-^(a;,z). 

The above observation in conjunction with Lemma 5 in [42] is sufficient 
to prove the equivalence of the sphere bundles C'j j ^ and Cjj^. But we need 
to prove a more general equivalence, involving all the sphere bundles C'j 
simultaneously, for < j < £. 

However, note that the proof of Lemma 5 in |42| proceeds by induction 
on the skeleton of the CW-complex of the base of the bundle. After choos- 
ing a sufficiently fine triangulation of the set Fj ^ = Fj ^ compatible with 
the closed subsets Fi ■ = Fi ■, the same proof extends without difficulty 
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to this slightly more general situation to give a fiber preserving homeomor- 
phism, (p : — > C^, which restricts to an isomorphism of sphere bundles, 
(ki,] '■ ~^ ^/j,z' each / C [m] and < j < ^. 

We have the following maps. 

(4.19) 

Ext(^x,R') Ext(Sx,R') Ext(Cx,R') ^ C'^ 

Ext(^z,R') Ext(Bz,R') Ext((:7z,R') ^ ^ C'^ 

The map « is the inclusion map, and r is a retraction shown to exist by 
Proposition 14.111 

Since all the maps 02, i, r have been shown to be homotopy equivalences, 
by Propositions 14.61 14.31 and 14. IH their composition is also a homotopy 
equivalence. 

Moreover, for each / C [m], the maps in the above diagram restrict 
properly to give a corresponding diagram: 
(4.20) 

Ext(^/,x,R') Ext(B,,,,R') Ext(C7,x,R') C;,, 

Ext(A/,z,RO Ext(S/,2,R') Ext(C/,z,R') ^ C\^^ 

Now let -Dx,z = C'x) ^^"^ /x = (j)2 ° "i ° f and = (j)2 o i o r o (p. Finally, 
for each / C [m], let -D7,x,z = C|x ^'^d the maps fi^:x.-,fi,x the restrictions 
of /x and /z respectively to -D/,x,z- The collection of sets Di^^^ and the 
maps //,x) //.z clearly satisfy the conditions of the lemma. This completes 
the proof of the lemma. □ 

Remark 4.25. Note that if Ri is a real closed sub-field of R, then Lemma 
14.241 continues to hold after we substitute "T C Ri" and "for all x G RJ^" 
in place of "T C R'^" and "for all x S R'^" in the statement of the lemma. 
This is a consequence of the Tarski-Seidenberg transfer principle. 

With the same hypothesis as in Lemma 14.241 we also have, 

Lemma 4.26. There exists a finite set T CK^ with 

#r < (2™£A;d)'^(™'=) 

such that for every x G R'^, there exists z £ T, for each I C [m], a semi- 
algebraic set -E'/,x,z defined over R" , where R" = R(e,e, ^) (see \4.lO\ for the 
definition of e and 5), and S-maps gi^x-,9i,x os shown in the diagram below 
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such that 5'7,xi5/,z are both stable homotopy equivalences. 
(4.21) 

Ext(f| ^i,x, R") Ext(f| Ai^,, R") 

For each / C J C [m], Sj,x,z C -E/.x.z and the maps 5j,x)5j,z are re- 
strictions of of gi^^^gi^^. 

Proof. Let 1 ^ e > be an infinitesimal. For 1 < -i < m, we define 

Q^ = Qi+e{Yi + --- + Yl), 

ii = {(y,x) I |y| = 1 A Q,(y,x) <0)}. 

Note that the set Aj^x is homotopy equivalent to Ext(f^ jdj^x, R(£)) 

is/ iG/ 

for each / C [m] and x S R*^ . Applying Lemma 14.241 (see Remark I4.25P 
to the family Q = {— Qi, . . . , —Q-m}^ we have that there exists a finite set 
T <ZR^ with 

#r < (2™£A;d)^('"*^) 

such that for every x S R*^, there exists z G T such that for each / C [m], 
the following diagram 

-D/,x,z 

Ext(|Ji,,x,R") Ext(|Jii,,,R") 

where for each x E R'"' we denote 

4x = {(y,x) I |y| = l A -Q,(y,x) < 0)}, 

//,X5 //,z are homotopy equivalences. 

Note that for each x G R'^, the set Ext(|^ ^i,x) R") is a deformation re- 

tract of the complement of Ext([^ ^j,x) R") and hence is Spanier- Whitehead 

dual to Ext(|^ Jij^x) R")- The lemma now follows by taking the Spanier- 

i&I 

Whitehead dual of diagram (j4.22p above for each / C [m]. □ 

Proof of Theorem I4.22[ Follows directly from Lemma [i?2il □ 

Proof of Theorem 14.231 Follows directly from Lemma [06l □ 
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We now prove a homogenous version of Theorem 14.11 

Theorem 4.27. Let R be a real closed field and let 

Q = {Qi,...,Qm}cR[Yo,...,Y,,Xi,...,Xk], 

where each Qi is homogeneous of degree 2 in the variables Yq, . . . ,Y£, and of 
degree at most d in Xi , . . . , Xk ■ 

Let vr : X R'^ — > R*^ 6e the projection on the last k co-ordinates. 
Then, for any Q-closed semi- algebraic set S C x R'^, the number of stable 
homotopy types amongst the fibers S'x is bounded by 

Proof. We first replace the family Q by the family, 

Q' = {Qu---,Q2m} = {Q,-Q I Qg Q}. 

Note that the cardinality of Q' is 2m. Let 

Ai = {{y,^) I |y| = 1 A Qi(y,x) <0)}. 
It follows from Lemma 14.261 that there exists a set T C R'^ with 

#r < (2'"£A;d)^('"*^) 

such that for every / C [2m] and x G R'^, there exists z G T and a semi- 
algebraic set -E/^x.z defined over R" = R{e,e,6) and S-maps ^/^xi^/.z as 
shown in the diagram below such that gi^x,gi,z are both stable homotopy 
equivalences. 



(4.23) 




Ext(f| R") Ext(f| ^i,,, R") 

Now notice that each Q-closed set S is a union of sets of the form |^ Ai 
with / C [2m] . Let 

s= u n^^- 

Moreover, the intersection of any sub-collection of sets of the kind. Hie/ ^« 
with / C [2m], is also a set of the same kind. More precisely, for any S' C S 
there exists G 21^™'] such that 
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We are not able to show directly a stable homotopy equivalence between 
5x and Sz- Instead, we note that the S-maps (7/,x and gi^z induce S-maps 
(cf. Definition EH]) 

5x : hocolim({Ext(p| ^i,x,R") U G S}) — > hocolim({£;7,x,z U G S}) 

Qz : liocolim({Ext(p| Aj,z,R") | / G S}) — > hocolim({£;/,x,z U G S}) 

which are stable homotopy equivalences by Lemma 12.331 since each gi^y^ 
and gi^z is a stable homotopy equivalence. 

Since hocolim({j^ ^j^x I I G S}) (resp. hocolim({f^ Aj^z | / G S})) is 

homotopy equivalent by Lemma [2.321 to ^i,x (resp. [^^j,z), it 

follows (see Remark [2.ip that 'S'x = ^i,x is stable homotopy equivalent 

to 5z; = Ai^z. This proves the theorem. □ 

/GSiG/ 

5.2. Inhomogeneous case. We are now in a position to prove Theo- 



rem 



Proof of Theorem 14.11 Let cphe a 'P-closed formula defining the V- 
closed semi-algebraic set S C R^^^'. Let 1 » e > be an infinitesimal, and 
let 



\i=i i=i / 



Let V = V L) {Po}, and let (j) be the "P-closed formula defined by 

<^ = 0A{Po<O}, 

defining the P-closed semi-algebraic set St C R(e)^^^. Note that the set 5;, 
is bounded. 

It follows from the local conical structure of semi-algebraic sets at in- 
finity |27j that the semi-algebraic set Sf, has the same homotopy type as 
Ext(5,R(e)). 

Considering each Pi as a polynomial in the variables Yi,...,Yi with 
coefficients in R[Xi, . . . ,Xk], and let P^^ denote the homogenization of Pi. 
Thus the polynomials Pj^ € K[Yq, . . . ,Y£, Xi, . . . , X^] and are homogeneous 
of degree 2 in the variables Yq,. . . ,Y£. 

Let C X R(e)'= be the semi-algebraic set defined by the P'^-closed 
formula (f)^ (replacing Pi by Pj^ in cf)). It is clear that S*^ is a union of two 
disjoint, closed and bounded semi-algebraic sets each homeomorphic to Sb, 
which has the same homotopy type as Ext(5, R(e)). 
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The theorem is now proven by applying Theorem 14.271 to the family 
and the semi-algebraic set S^. Note that two fibers S-x_ and Sy are 
stable homotopy equivalent if and only if Ext(5'x, R(e)) and Ext(S'y, R(e)) 
are stable homotopy equivalent (see Remark 12. ip . □ 

6. Metric upper bounds 

In |22| certain metric upper bounds related to homotopy types were 
proven as applications of the main result. Similar results hold in the qua- 
dratic case, except now the bounds have a better dependence on I. We state 
these results without proof. 

We first recall the following results from |22j . Let y C R^ be a V-semi- 
algebraic set, where V C Z[yi, . . . , Y^]. Suppose for each P gV, deg(P) < d, 
and the maximum of the absolute values of coefficients in P is less than some 
constant M, < M G Z. 

Theorem 4.28. There exists a constant c > 0, such that for any ri > 
r2 > M'^"^ we have 

(1) yni?^(0,ri) and V r\ Bi{d,r2) are homotopy equivalent, and 

(2) V\Bi{d,ri) and V \ Bi{Q,r2) are homotopy equivalent. 

In the special case of quadratic polynomials we get the following im- 
provement of Theorem 14.281 

Theorem 4.29. Let R be a real closed field. Let V <ZK^ he a V-semi- 
algebraic set, where 

V = {Pi,...,Pm}cR[Y^,...,Ye], 

with deg(Pj) < 2, 1 < i < m and the maximum of the absolute values of 
coefficients in V is less than some constant M , < M € Z. 

There exists a constant c > 0, such that for any r\ > r2 > M^""" we 
have, 

(1) V r\Bi[{),ri) andVr\Bi{Q^r2) are stable homotopy equivalent, and 

(2) V \ Bi{0,ri) and V \ Bi{0,r2) are stable homotopy equivalent. 
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Algorithms and Their Implementation 

1. Computing the Betti Numbers of Arrangements 

In this chapter, we consider arrangements of compact objects in M'^ which 
are simply connected. This impHes, in particular, that their first Betti num- 
ber is zero. We describe an algorithm for computing the zero-th and the 
first Betti number of such an arrangement, along with its implementation 
[15j . For the implementation, we restrict our attention to arrangements in 

and take for our objects the simplest possible semi-algebraic sets in 
which are topologically non-trivial - namely, each object is an ellipsoid de- 
fined by a single quadratic equation. Ellipsoids are simply connected, but 
with non- vanishing second co- homology groups. We also allow solid ellip- 
soids defined by a single quadratic inequality. Computing the Betti numbers 
of an arrangement of ellipsoids in is already a challenging computational 
problem in practice and to our knowledge no existing software can effectively 
deal with this case. Note that arrangements of ellipsoids are topologically 
quite different from arrangements of balls. For instance, the union of two 
ellipsoids can have non-zero first Betti number, unlike in the case of balls. 

1.1. Outline of the Method. The following corollary follows imme- 
diately from Proposition 12.201 

Corollary 5.1. Let be S = UilLi S'j C M.'^ such that Si,...,Sm are 
compact semi- algebraic sets with 

(1) R^{Si) =Q, and 

(2) il\Si) = 0, 1 < i < m. 

Let the homomorphisms Sq and 5i in the following sequence be defined as in 
Chapter \2.2\ (identifying H'^(K) with the Q-vector space of locally constant 
functions on a simplicial complex K). 

e. hO(5,)^ e.<, H0(5. n Sj)^ e,<^.<, r^{S, n s, n Se). 

Then, 

bo{S) = dim(Ker(<5o)), 

bi{S) = dim(Ker(5i)) -dim(Im((5o)). 

The importance of Corollary 15.11 lies in the following observation. Given 
an arrangement, . . . ,Sm}, of m simply connected objects in M^, sup- 
pose we are able to identify the connected components of all pairwise and 

65 



66 



5. ALGORITHMS AND THEIR IMPLEMENTATION 



triple-wise intersections of these objects and their incidences (that is, which 
connected component of Si Ci Sj D Sg is contained in which connected compo- 
nent of Si n Sj). Then this information is sufficient to compute the zero-th 
and the first Betti number of the arrangement. We only have to look at 
the objects of the arrangement at most three at a time. Thus, the cost of 
computing the connected components and incidences is O(m^). This is to 
be compared with having to compute a global triangulation of the whole 
arrangement using cylindrical algebraic decomposition which would have 
entailed a cost of 0{m ). 

Recall that a cylindrical decomposition (see Chapter II. 4p adapted to a 
finite set V of polynomials in M[Xi, . . . , Xk\ produces a graph where the ver- 
tices correspond to cells in Sk and edges correspond to adjacencies. More- 
over, each cell in Sk is P-invariant and we know the sign for each P in 
V on each such cell. Hence, given an arrangement, {S'l, . . . , 5m}, of m 
semi-algebraic sets in M*^, we are able to identify the connected compo- 
nents of all pairwise and triple-wise intersections of these objects and their 
incidences by computing a cylindrical decomposition adapted to the fami- 
lies Vij/, 1 < i < j < £ < m, where Vij^i is the set of polynomials used 
in the definition of Si,Sj, and Se and by performing a graph transversal 
algorithm on the graph described above. 

To sum up, we now formally describe our algorithm for computing the 
zero-th and the first Betti numbers of an arrangement of m simply connected 
compact objects in M^. 

Algorithm 5.2 (Computing the zero-th and the first Betti number). 
Input: compact sets Si C M^', I < i < m, with bo{Si) = 1 and 6i(5j) = 0. 
Output: bQ^S) and bi{S). 
Procedure: 

• For each triple 1 < i < j < ^ < m, do the following: 
Compute a cylindrical decomposition adapted to the set {Si, Sj, Sg}. 
Identify the connected components of all pairwise and triple-wise 
intersections and their incidences. 

• Compute the matrices A and B corresponding to the sequence of 
homomor phisms : 

H0(5,)^ e,<, B^s, n s,)^ e^<J<i hO(5, n s, n Se). 

• Compute 

60(5*) = do — r:k{A), and 
bi{S) = di-Tk{B) -ik{A), 

where is the dimension of 0i<j<m H'^('S'i)) di is the dimension 
of 

®i<i<j<m^^i'^i^ '^j)' the rank of a matrix is denoted by 

rk(.y. 
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1.2. The Implementation. The algorithm has been prototypicahy 
implemented using QEPCAD B (Version 1.27) |66j and Magma |60j for 
compact sets C M^. We use the package QEPCAD B for computing the 
cylindrical decompositions, in Step 1 of Algorithm 15. 2[ There are several 
other packages available for computing cylindrical decompositions, for in- 
stance REDLOG |67j . The main reason for using QEPCAD B is that it 
provides some important information regarding cell adjacency, that is not 
provided by the other systems. 



— (l)pl(-,-) 


— {l,l)pl(+) 


(l,l.l)pl(+) 


F 


— (2)pl(0,-) 


— {2,l)pl(+) 


— (2,l,l)pl(+) 


F 




(2,2)pl(0) 


(2,2,l)pl(+) 


F 






— (2,2,2)pl(0) 


T 






— (2,2,3)pl(+) 


F 




— (2,3)pl(+) 


(2,3,l)pl(+) 


F 


(3)pl(+,-) 


— (3,l)pl(+) 


(3,l,l)pl(+) 


F 




— (3,2)pl(0) 


(3,2,l)pl(+) 


F 






— (3,2,2)pl(0) 


T 






— (3,2,3)pl(+) 


F 




— (3,3)pl(-) 


(3,3,l)pl(+) 


F 






— (3,3,2)p2(0) 


T 






— (3,3,3)pl(-) 


F 






— (3,3,4)p2(0) 


T 






— (3,3,5)pl(+) 


F 




— (3,4)pl(0) 


— (3,4,l)pl(+) 


F 






— (3,4,2)pl(0) 


T 






— (3,4,3)pl(+) 


F 




— (3,5)pl(+) 


(3,5,l)pl(+) 


F 


— C4)pl(+,0) 


— (4,l)pl(+) 


(4,l,l)pl(+) 


F 




— (4,2)pl(0) 


— (4,2,l)pl(+) 


F 






— (4,2,2)pl(0) 


T 






— (4,2,3)pl(+) 


F 




— (4,3)pl(+) 


(4,3,l)pl(+) 


F 


— (5)pl(+,+) 


— {5,l)pl(+) 


— (5,l,l)pl(+) 


F 



Figure 1. Output of a cylindrical decomposition using QEPCAD B 

In Figure [H which shows the QEPCAD B output for a cylindrical de- 
composition adapted to the unit sphere, the first (resp., second and third) 
column corresponds to the cylindrical decomposition of the line (resp. plane 
and M'^). Note that the signs accompanying the cells give the signs of pro- 
jection factors computed by QEPCAD B and the letter "T" and "F" corre- 
sponds to true and false value of the cells, i.e., depending upon whether our 
input formula is true or false on this cell. 

Even though QEPCAD B does not provide full information regarding 
cell adjacencies in dimension three, we are still able to deduce all the needed 
cell adjacencies as described in Chapter 11.4.31 rnaking use of the fact that 
input polynomials are quadratic. 

We use Magma for post-processing of the information output by QEPCAD B, 
in Steps 2 and 3 of the algorithm. Note that all computations performed 
are exact with no possibility of numerical errors. 

To illustrate our implementation, we consider four examples where the 
ellipsoids 

5i = {(x,y,z) I Pi(x,y,z) =0}, 
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1 < i < 27, are defined by the following list of polynomials (see Table [I]) 
We denote by A and B the matrices of the homomorphisms 6i and 82 with 

Table 1. Input polynomials defining the different arrangements 







8/9Xf + 1/64X| + 1/6X1 - 1 




P2 




l/64Xi2 + 8/9X| + 8/9X| - 1 




-rs 




8/9Xf + 8/9X| + 1/64X| - 1 








8/9(Xi - 4)2 + 1/64(X2 - 4)2 + 1/6X| 


-P5 




l/64(Xi - 4)2 + 8/9(^2 - 4)2 + 8/9X| 


-re 




8/9(Xi - 4)2 + 8/9(X2 - 4)2 + 1/64X| 






[Xi - ly + (X2 - 2)^ + X| - 3 








\^2 1 1 /n \^2 i o \^2 i 
5A{ + i/yA2 + 2A| — i 








l/9Xi2 + 5X| + 5X32 - 1 




-rio 




5Xi2 + 5X| + l/9X|-l 




D 

-Til 




5(Xi - 1)2 + 1/9(^2 - 1)2 + 2X32 


- 1 


D 

-ri2 




l/9(Xi - 1)2 + 5(^2 - 1)2 + 5X32 


- 1 


-ri3 




5(Xi - 1)2 + 5(^2 - 1)2 + 1/9X32 


- 1 


P, . 
-ri4 




5{Xi + 1)2 + 1/9(^2 - 1)2 + 2X32 


- 1 


As 




l/9(Xi + 1)2 + 5(^2 - 1)2 + 5X32 


- 1 


Ae 




5(Xi + 1)2 + 5(^2 - 1)2 + 1/9X32 


- 1 


At 




5{Xi - 1)2 + 1/9(^2 + 1)2 + 2X32 


- 1 


As 




l/9(Xi - 1)2 + 5(^2 + 1)2 + 5X32 


- 1 


A9 




5(Xi - 1)2 + 5(X2 + 1)2 + 1/9X32 


- 1 


Ao 




5(Xi + 1)2 + 1/9(X2 + 1)2 + 2X32 


- 1 


Ai 




l/9(Xi + 1)2 + 5(X2 + 1)2 + 5X32 


- 1 


A2 




5(Xi + 1)2 + 5(X2 + 1)2 + 1/9X32 


- 1 


As 




6(Xi - 1/2)2 ^ + i/Qxl - 1 




A4 




4X^2 + 4(X2 - 1/2)2 + i/Qxl - 1 




As 




5(Xi + 2)2 + 5X1 + 1/6X32 - 1 




Ae 




l/6(Xi + 2)2 + 5(X2 - 2)2 + 5X32 


- 1 


A7 




5(Xi + 2)2 + 1/6(X2 - 2)2 + 5X32 


- 1 



respect to the obvious basis. The columns (resp., the rows) of the matrix 
A are labeled by (resp., j), while the columns (resp., the rows) of the 
matrix B are labeled by ef • (resp., ef .,), where corresponds to Si, ef- 
corresponds to the p-th connected component of SiCiSj and e^j ^ corresponds 
to the p-th connected component of Si H Sj n Si. 

Remark 5.3. In the examples described below, we have modified the 
matrix A as follows. Since we know that each input set Si has exactly one 
connected component, we can simplify the computation. We only need to 
check whether or not the intersection Si D Sj is empty. Therefore, we have 
exactly one row for each intersection instead of one row for each connected 
component of each intersection SiCiSj, and this reduces the size of the matrix 
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A without changing its rank. For the matrix B we delete all rows containing 
only zeros which correspond to empty triple intersections Si n Sj H Si. 




Example 5.4 (Three ellipsoids). Let S be the union of the first three 
ellipsoids, i.e., S = IJj=i (see Figure [2]). Then 



A 



ei 62 63 



-1 1 


ei,2 


-1 1 


ei,3 


{ -1 I ) 


62,3 



B 



1 


pi 

^1,2 


^1,2 


ei,3 


62,3 


\ 




1 





-1 


1 






1 





-1 


1 






1 





-1 


1 






1 





-1 


1 









1 


-1 


1 









1 


-1 


1 









1 


-1 


1 




\ 





1 


-1 


1 
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In this case, 



6o(5) = do - rk(A) = 3-2 = 1 

bi{S) = (ii - rk(B) - rk(A) = (4 - 2) - 2 = 



Example 5.5 (Six ellipsoids). Let the set S be the union of the first six 
ellipsoids Sj, 1 < i < 6, i.e., S = Ui=i Si (see Figure [3]). Then 
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bo{S) 
bi{S) 



= do- rk{A) =6-5 = 1 

= di- rk{B) - rk(^) = (12 - 4) - 5 = 3 
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Figure 4. Seven ellipsoids 



Example 5.6 (Seven ellipsoids). Let the set S be the union of the first 
seven ellipsoids S",, 1 < i < 7, i.e., S = (JLi (see Figured]). Then 
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In this case, 
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6o(5) = do - rk(A) = 7-6 = 1 

hx{S) = di - rk(5) - rk(^) = (16 - 7) - 6 = 3 




Figure 5. Twenty ellipsoids 



Example 5.7 (20 ellipsoids). Let the set 5 be the union of the last 20 
ellipsoids Si, 8 < i < 27, i.e., S = IjfJ^gSi (see Figure [5]). Thus, we get a 
190 X 20-matrix A of rank equal to 19, a 190 x 107-matrix B of rank equal 
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to 55, and the dimension of 0j<j H^{Si n Sj) is equal to 107. In this case, 
6o(5) = do - rk(^) = 20 - 19 = 1 
bi (S) = di- rk{B) - rk(A) = 107 - 55 - 19 = 33 

2. Computing the Real Intersection of Quadratic Surfaces 

In this chapter, we consider the problem of computing the real inter- 
section of three quadratic surfaces, or quadrics, defined by the quadratic 
polynomials Pi, P2 and P3 in M^. We describe an algorithm for computing 
the isolated points and a linear graph embedded into (if the real inter- 
section form a curve) representing the real intersection of the three quadrics 
defined by the three polynomials Pi, along with its implementation |52| . 
For the implementation, we restrict our attention to quadrics with defining 
equation having rational coefficients. 

Before outlining our method, we define the silhouette curve and cut 
curve, which can be interpreted in our setting as the projection of one 
quadric and the projection of intersection curve of two quadrics into the 
Xi — X2-plane, respectively. 

Definition 5.8. Let P, Q £ R[Xi, X2, X3]. The algebraic curves with 
defining polynomials 

Sil(P) := Res(P, dP/dXs), cut(P, Q) := Res(P, Q) 

are called silhouette curve and cut curve respectively. 

Another geometric interpretation of the silhouette curve in our set- 
ting is the following. The silhouette curve defined by Sil(Pi) contains all 
points (x, y) such that the polynomial Pi(x, y,X3) has exactly one root z 
of multiplicity 2. 

2.1. Outline of the Method. The basic idea of computing the in- 
tersection of three quadrics is based on the cylindrical decomposition (see 
Chapter II. 4p . As the algorithm of Schomer and Wolpert |69|i I83| , our ap- 
proach can be summarized by several phases: 

preparation, projection, planar arrangement analysis and lifting 
phase. 

But our analysis of the planar arrangement and lifting phase differs from 
the methods presented in |69|, I83j . 

First, we project one input quadric and the resulting space intersection 
curves of the pairwise intersections by computing (univariate) resultants 
onto the plane assuming that we have a "good" coordinate system by us- 
ing the Brown- McCallum projection operation (see [301, I29j ). The Brown- 
McCallum projection operation produces, based on the current literature, 
the smallest projection set in our setting. Then we analyze the planar ar- 
rangement of curves before we lift our solution into space (if possible). In 
other words, we compute the defining polynomial Sil(Pi) of the silhouette 
curve of the input quadric Pi, and the defining polynomial cut(Pi, Pj) of the 
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corresponding cut curves (see Definition I5.8p . Then we identify tfie com- 
mon factor G of cut(Pi, P2) and cut(Pi, P3) and the corresponding gcd-free 
parts Hi = cut{Pi, Pi)/G. 

While Schomer and Wolpert |69j use resultants for computing the candi- 
dates of the Xi and X2-coordinates and analyze the resulting grid (see |69j 
for more details) afterwards, our planar analysis is based on the TOP algo- 
rithm (see Algorithm 12. 12p . To be more precise, we use the TOP-algorithm 
in order to obtain the topology of the curve defined by the common factor G 
(including some relevant points on the curve). In addition, we generalize the 
idea of |45| of using subresultants to two planar curves. We perform a linear 
change of coordinates (if it is needed) in order to have two planar curves of 
the arrangement in generic position (see Definition l2.10p . Then we compute 
the Xi-coordinates of the (planar) intersection points of two curves using 
resultants as well as the X2-coordinates of those points that can be described 
rationally in terms of Xi via subresultant computations. 

Furthermore, when the intersection points form a curve, the set of solu- 
tions is described topologically via a linear graph embedded in M'^. The com- 
puted graph provides all the information for tracing this curve numerically 
since we know exactly how to proceed when we are close to a complicated 
point. Nevertheless, all computed points lie in the real intersection set of 
the three quadric surfaces defined by the three polynomials Pi. 

Next we describe all necessary steps, but we omit the details on the 
change of coordinates to which we refer to several times. For more details 
on the change of coordinates see |45j . 

2.2. Details on the Preparation Phase. Before starting the real 
computation, we test the input quadrics for degeneracy and if they behave 
well under resultant computations. We want to make sure that all quadrics 
are of degree equal to two, Xs-regular, square-free and pairwise do not have 
a common factor of degree equal to one. The absence of the latter two 
conditions can easily be detected and solved as it simplifies the considered 
problem. For example, one quadric describes a single plane if it is not square- 
free, whereas two quadrics define three planes, with one of them in common, 
if they have a common factor of degree equal to one. Therefore, we omit 
the details on these cases, though we can detect and solve them easily. In 
the case that one of the other conditions is violated, we make a change of 
coordinates and start the computation again. 

Finally, from now on we use the following assumption for the input 
polynomials Pi. 

Assumption 5.9. The trivariate polynomials Pi, P2 and P3 with co- 
efficients in M are all of degree equal to two, square-free, Xa-regular and 
pairwise do not have a common factor of degree equal to one. 

2.3. Details on the Projection Phase. After a suitable preparation 
of our input quadrics, we assume throughout this and the following sections 
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that the quadratic input polynomials have the properties of Assumption l5.9l 
It is worthwhile to mention that Assumption 15.91 is necessary in order to in- 
terpret correctly the projection onto the X1-X2 plane via resultant computa- 
tion. Our projection method is based on the so-called restricted equational 
version of the Brown-McCallum projection operation (see |30j . |29j ) where 
we use the polynomial Pi as the pivot constraint. The Brown-McCallum 
projection operation produces, based on the current literature, the smallest 
projection set V and consists of the following polynomials in our setting, 

V = {Sil(Pi),CUt(Pi,P2),CUt(Pi,P3)}. 

As in the beginning of our computation, we need to test the polynomials 
contained in the set V for degeneracy in order to interpret correctly the 
following resultant computations. Thus, we simplify the set V further and 
we obtain the set V containing the following polynomials, 

P = {Sil(Pi),if2,if3,G}, 

such that Hi = cut(Pi, where G is the greatest common divisor of 
cut(Pi,P2) and cut(Pi,P3). Moreover, we decompose the polynomial G 

further. We write G = G ■ Sil(Pi) where G (resp., Sil(-Pi)) is the gcd- 
free part (resp., greatest common divisor) of G and Sil(Pi). Note, that the 

decomposition of the polynomial G into Sil(Pi) and G will be very useful for 
the lifting phase (see Chapter l2.5p . Finally, we can summarize the projection 
phase as follows. 

Algorithm 5.10 (Projection). 
Input: three polynomials Pi, P2 and P3 in M[Xi,X2,X3] with the proper- 
ties of Assumption 15. 9i 
Output: V = {Sil(Pi),i/2,^3,G} 

such that Hi is the square-free part of cut(Pi, Pi) with respect to G, 
where G = gcd(cut(Pi, P2), cut(Pi, P3). Moreover, we decompose the poly- 
nomial G into G = G ■ Sil(Pi) where G (resp., Sil(Pi)) is the gcd-free part 
(resp., greatest common factor) of G and Sil(Pi). 

2.4. Details on the Analysis of the Planar Arrangement. In this 
section, we describe how we analyze the planar arrangement. We assume 
from now on that the set V computed before is of the following form: 

V = {Sil{Pi),H2,Hi,G}, 

such that Hi = cut(Pi, Pj)/G, where G is the square-free part of the com- 
mon factor of cut(Pi,P2) and cut(Pi,P3)). Moreover, we decompose the 

polynomial G further. We write G = G ■ Sil(Pi) where G (resp., Sil(Pi)) is 
the gcd-free part (resp., greatest common factor) of G and Sil(Pi). 

The problem, which might occur, is that the planar curves might not 
be in generic position which would ensure that we can use subresultants 
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in order to compute the critical points (including intersection points with 
another curve) of the planar curves in our arrangement. In this case, we start 
the computation again after a change of coordinates if the planar curves are 
not in generic position. Hence, we assume throughout this and the following 
sections that the set V has the following properties. 

Assumption 5.11. Let V = {S\\{Pi),H2,H^,G] as computed in Algo- 
rithm 15.101 such that all polynomials are X2-regular. The polynomials H2 
and Hj, as well as G are in generic position. Moreover, G is in generic 
position with respect to Sil(Pi). 

By using the Brown-McCallum projection operation for eliminating the 
variable X3, it follows that the (possible) intersection points of all three 
quadrics lie on the cut curves defined by cut(Pi,P2) and cut(Pi,P3), i.e., 
on the intersection of Zer(i?2,K^) and Z(iT:{H^,R^), or on Zer(G,M^). In 
addition, we need to identify the common points of those curves with the 
silhouette curve Zer(Sil(Pi), M^) since the number and type of points above 
a point on Zer(Sil(Pi), M^) might be different than for points which do not 
lie on Zer(Sil(Pi), M^). But observe that the curve Zer(Sil(Pi), M^) contains 
all points (x, y) such that -Pi(x, y, X3) has exactly one root z of multiplicity 
2. To sum up, we need to compute the following: 

(1) the intersection points of Zer(i/2) 1^^) and Zev{H^, M?) and whether 
or not they lie on the curve Zer(Sil(Pi), M^), and 

(2) the topology of Zeic{G, M^) including the common points with Zer(Sil(Pi), M^) 
which could be finitely or infinitely many. 

By decomposing the polynomial G = G ■ Sil(Pi) we simplify the second 
problem further since we just need to compute the following: 

1) the topology of Zer(G',M^^-including the common points with Sil(Pi), 

2) the topology of Zer(Sil(Pi)) including the common points with 

Zer((5,M2). 

It is worthwhile to mention that we can not decide without further compu- 
tation whether or not a planar point can be lifted to a solution of all three 
quadrics. This comes from the fact that two different (space) points in 
might get projected to the same (planar) point. Nevertheless, this problem 
can be solved easily as we will see in Chapter 12.51 We summarize the above 
discussion in the following algorithm. 

Algorithm 5.12 (Planar Arrangement Analysis). 
Input: the set of polynomials 

V = {Sil{Pi),H2,H3,G}, 

with the properties of Assumption 15.111 
Output: 

• the common points of Zer(iJ2,R^) and Zer(iJ3,M^), 

• the topology of the curve Zer(G,M^), described by 
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— the real roots xi,. . . ofRes{G, 80/8X2), Res(Sil(Pi), 5Sil(Pi)/5X2) 
and _ 

Res(G, Sil(Pi)). We denote by xq = —00, x^+i = 00. 

— The number nii of roots of G(x, X2) in M when x varies on 
(xj,Xj_|_i). We denote this root by Xj^i, . . . ,Xj^m-. 

— The number rii of roots of G{xi,X2) in M. We denote these 
roots by yi,i, . . . ,yj,„,. 

— A number q < rii such that if (xi,Zj) is the unique critical 
point of the projection of Zer(G, C^) on the Xi-axis or an 
intersection point of Zer(G,M^) and Zer(Sil(Pi), M^) above Xj, 

= yi,Ci ■ 

Procedure: 

- Compute the common points of Zer(f/'2) I^^) and Zer(if3,M^) using Algo- 
rithm [2l2] (TOP) as a black-box. 

- Compute the topology of Zer(G, M^) including the common points with 
Zer(Sil(Pi),M2) by using Algorithm [232] (TOP) as a black-box. 

2.5. Details on the Lifting Phase. 

2.5.1. Lifting of Single Points. We recall some well-known facts about 
the real roots of a quadratic polynomial in one variable. Note that we know 
a priori what case we do have to consider in our setting. For example, a 
candidate (x, y) G Zer(G, M^) which also lie on Zer(Sil(Pi), M^) corresponds 
to the case D = (see Proposition I5.13p . i.e., the polynomial Pi{x,y, X^) 
has exactly one real root. 

Proposition 5.13. Let P = aX"^ + bX + c with a, 6, c G M and let D = 
6^ — 4ac. Then we get the following cases. 

(1) If D = 0, then the polynomial P has exactly one solution x = — 

(2) If D > 0, then the polynomial P has two real solution xi and X2. 

In this case, xi = ^ (^—b — ^/D^ and X2 = ^ (^—b + \fD^ . 

(3) If D < 0, then the polynomial P has only two complex conjugated 
roots. 

By using the information computed by Algorithm 15 . 1 2 1 we can now easily 
determine the solutions zi,...,Zj, i < 2, of the polynomial Pi{x,y, X^) 
where (x, y) is a possible candidate in the plane. 

2.5.2. Lifting of a Curve. Our approach for lifting a curve is similar 
to lifting a single point as described in the chapter before. By computing 
some extra points on Zer(Pi,M^) as described in the previous section, we 
can determine easily the adjacency of the (possible) space curve which is 
induced by the plane curve Zer(G, M^). 

Assume that we computed the topology of Zer(G, M^) as described in 
Algorithm 15. 12[ First, we lift all points (if possible) onto Zer(Pi,M^). Note 
that we can easily determine the missing adjacencies as described in Chap- 
ter II. 4. 3^ since there are only one or two points above. Then we just need 
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to test whether or not our candidates he on Zer(i-2) 1^'^) and Zer(P3,]R^) 
as well. It is worthwhile to mention that not all components of Zer(G, M^) 
might get lifted even though they can be lifted to a solution on Zer(Pi, R^). 



2.6. The Implementation. The algorithm has been prototypically 
implemented in the Computer Algebra System Maple (version 9.5) [61] and 
it follows the approach outlined closely. It starts always with three quadratic 
polynomials Pi, P2 and P3 in Q[Xi, X2, X3] and, due to efficiency reasons, it 
performs most of the computations by using floating point arithmetic. The 
latter one comes from the fact that we extended Laureano Gonzalez- Vega 
and loana Necula's TOP algorithm code ( [45] ). Hence, the only computa- 
tions that are performed symbolically are: 

(1) the computation of the projection set V = {Sil(Pi), H2,H3, G}. 

(2) the computations of the different signed subresultant sequences and 
their coefficients for the projection set V. 

(3) the computation of the square-free part of the resultant of two 
polynomials Pi and P2 in Q[Xi,X2] and its decomposition with 
respect to the signed subresultant coefficients. 

The remaining computations consist in solving numerically different poly- 
nomial equations (without multiple roots) or evaluating at these roots some 
of the polynomials symbolically computed. Initially the chosen precision is 
15 digits, but one can choose any other starting precision ti. As in the im- 
plementation of the TOP-algorithm, we choose a threshold e that depends 
on the chosen precision in order to decide whether or not a polynomial is 
zero at a given point. 

Once the planar arrangement V is computed, we analyze the size t2 
of the input polynomials Pi and the set V. Afterwards, we update the 
precision to t digits, where t = max{ti,t2 + 10,15}. Furthermore, the 
Maple function f solve is used to solve the square-free univariate polynomial 
equations before mentioned. If f solve does not return the correct number 
of roots (which are known in advance) or some numerical evaluation returns 
some non guaranteed value, the precision is increased by 10 digits and those 
computations are performed again. Moreover, we output the coordinates of 
the isolated points and a three dimensional linear graph if the intersection 
points form a curve. 

We end this section by giving some examples, which illustrate our ap- 
proach. The experimentations were performed on a PowerPC G4 IGHz. 
The following example is taken from |69j . 
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Example 5.14 (two isolated points, [69]). Let be 

Pi = 7216Xf - 11022X1X2 - 12220X1X3 + 15624X| + I5I68X2X3 + 11186X| - 1000 
P2 = 4854Xf - 356OX1X2 + 4468X1X3 + 658Xi + 5040X| + 32X2X3 + 1914X2 + 

10244X| + 3242X3 - 536 
P3 = 8877Xf - 10488X1X2 + 9754X1X3 + I28OX1 + 16219X| - 16282X2 X3 - 8O8X2 + 

10152X| - III8X3 - 796 

Then the projection set V contains of 

Sil(Pi) = 10846519Xf- 7653903XiX2 - 2796500 + 29313252X1 

H2 = -56556109351696Xi + 61135807177688X2 - 6220192626724 + 

203315497528241XiX2 - 56404750618857X1^ - 55861103592035X| 
-910824371936818X2 Xi^ + 972629091137652XiX| - 659086885094112X| + 
533601 199106972Xf - 2885241224346328Xf X2 + 4223689039107028X1^X1 
-3571456229045952XiX| + 1026392565603269X^ + 1407622740362496X| 
F3 = 2872582087600Xi - 4005061111776X2 +69677228486124XiX2 

-23228971077672Xi^ - 49611754602456X| - 5464061993528X2 Xi^ 
-17976875889356XiX| + 40411859296976X| + 1462282618132Xf 
-926282674085672Xf X2 + 1733300718748310Xf X| - 1854003852157600XiX| + 
225439274765947Xf + 897407958763127X| - 66086625728 
G = 1 

Our computations end with a precision of 26 digits. The real intersection is 
computed in 0.572 seconds and consists of two isolated points, namely, 

/ -0.47111071472741316264056772 \ 
pi = -0.19897789206886601999604553 
\ 0.18592931583225857372754588 / 

and 




-0.16627634657169906116678201 \ 
0.10827914469994312737865267 
-0.011248383019525287650192532 / 



Moreover, Tableland Table [3] present a comparison between the com- 
puting times (in seconds) obtained by our approach and the prototypically 
and improved implementation of |69j^ using different numbers of decimal 
digits for the three input quadrics. Moreover, Table [2] contains the following 
additional information: 

• size of Input (resp., V) ~ number of decimal digits of the Input 
(resp., the projection set V. 

• Changes - number of linear changes of variables 

• Precision - used precision for obtaining the result 
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Table 2. Experimental results for Example 15.141 



Size of Input 


Size of V 


Changes 


Precision 


Time 


5 


16 


1 


26 


0.572 


8 


32 





42 


0.466 


12 


46 





56 


1.120 


15 


54 





64 


4.453 


20 


75 





85 


4.662 


23 


90 





100 


7.361 


28 


106 





116 


6.479 


32 


122 





132 


6.665 


36 


137 





147 


8.077 


40 


147 





157 


7.609 



Table 3. Experimental results of Schomer and Wolpert |69j 



Number 


of digits 


5 


10 


15 


20 


25 


30 


Running 


time 1 


18 


33 


56 


92 


126 


186 


Running 


time 2 


1.1 


2.7 


5.0 


7.8 


12.1 


16.1 



It is worthwhile to mention that we obtain similar running times for all 
our experiments. Additionally, the improvement of the running times do 
not only depend on the newer computer. 

Example 5.15 (closed curve). Let be 

Pi = {Xi-X2f + xl + xl-i 

P2 = {Xi-X2-lf + Xl + Xl-1 

Ps = 4X1 + 4X1-3 

Note, that the three quadrics are linearly independent and the projection 
set V contains of 

Sil(Pi) = - 2X1X2 + 2X| - 1 
H2 = 1 

Hs = -1-2X1 + 2X2 

sn(K^) = 1 

G = 1 - 2X1 + 2X2 

Then the real intersection of the three quadrics defined by Pi, P2 and P3 
consists of infinitely many points. Figure [6] shows the linear three dimen- 
sional graph computed by our implementation. The computations start and 
end with a precision of 15 digits and is computed in 0.101 seconds. For 
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Figure 6. The intersection of three hnearly independent quadrics 

representing the Unear graph we computed the following four points. The 
points 

= (-0.366025403784439,-0.866025403784439,0) 
P2 = (1.36602540378444,0.866025403784440,0) 

which correspond to the lift of the intersection points of the two plane 
curves Zer(Sil(Pi), M^) and Zer(G,M2), and 

(0.500000000000000, 0, -0.866025403784440), 
0.500000000000000, 0, 0.866025403784440) 

which are two sample points for the two curve segments between the critical 
points pi and p2 ■ 



Table 4. Experimental results for Example 15.151 



Size of Input 


Size of V 


Changes 


Precision 


Time 


0.101 


1 


1 





15 


0.185 


4 


8 





18 


0.257 


8 


15 





25 


0.196 


12 


20 





30 


0.307 


16 


30 





40 


0.323 


20 


38 





48 


0.498 


25 


47 


2 


57 


0.520 


28 


53 


2 


64 


0.591 


33 


62 


2 


82 


0.368 


36 


66 
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Example 5.16 (2 isolated points, G / 1). Let be 
Pi = 27X1^ + Q2Xl + 249X| - 10 

P2 = 88Xf + 45X| + 67X| - 66X1X2 - 25X1X3 + I2X2X3 - 24Xi + 2X2 + 29X3 
P3 = 88X^ + 45x1 + 67x1-66X1X2 + 25X1X3-12X2X3-24X1 + 2X2-29X3 

Note, that P3(Xi,X2,X3) = P2(Xi, X2, -X3). Then the projection set V 
contains of 

Sil(Pi) = 27Xf + 62X| - 10 

H2 = F3 = sn(A) = i 

G = -1763465 + 408332484X^ + 51939673X| + 10482900Xi - 2305740X2 
-123026916XiX| + 221120964X^X2 + 4764152X| + 17767644X| + 
14441004X1X2 - 250019406Xf + 16691919X1^ - 664779204X^X2 + 
564185724X^2^1 - 241015068XiX| 

Our computations end with a precision of 19 digits. The real intersection 
consists of two isolated points 

(0.06676451891748808143, 0.3991856119605212449, 0), 
(0.4954772252006942431, 0.2331952878577051550, 0) 

and is computed in 0.490 seconds. 



Table 5. Experimental results for Example 15.161 



Size of Input 


Size of V 


Changes 


Precision 


Time 


2 


9 





19 


0.490 


6 


22 





32 


0.355 


10 


37 





47 


2.374 


14 


46 





56 


4.939 


18 


67 





77 


5.018 


21 


82 





92 


6.362 


26 


98 





108 


6.515 


30 


113 





123 


7.109 


34 


129 





139 


7.694 


38 


138 





148 


9.671 


41 


158 





168 


9.056 



Example 5.17 (empty intersection). Let be 

Pi = X2 + Xi^ + 2X1X2 + 2X1X3 + X| + 2X2X3 + X| 
P2 = X| + 1 — X2 
P3 = 2x1 + 2-2X2 
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Then the projection set V contains of 
Sil(Pi) = X2 

H2 = i73 = sn(A) = i 

G = 1+Xf + Xl+4XfX2 + QXfxl + 4XiX'i- 4^X2 + 6X^ + 4^X1X2 + 2Xf 

Our computations start and end with precision of 15 digits. The real inter- 
section is empty and computed in 0.182 s. 



Table 6. Experimental results for Example 15.171 



Size of Input 


Size of V 


Changes 


Precision 


Time 


1 


1 





15 


0.182 


4 


15 





25 


0.191 


8 


30 





40 


0.187 


12 


39 





49 


0.274 


16 


59 





69 


1.025 


20 


74 





84 


0.978 


24 


92 


2 


121 


2.345 


28 


105 


1 


126 


1.863 


32 


122 


1 


142 


1.821 


36 


133 


1 


153 


2.090 


40 


152 


2 


182 


2.740 



Example 5.18 (a curve and an isolated point). Let be 

Pi = X2 + xf + 2x1X2 + 2x1X3 + xi + 2x2X3 + xi 
P2 = X3 — X2 + X1X2 + X2 + X2X3 
P3 = 2x1-2x2 + 2x1X2 + 2x1 + 2x2X3 

Note, that P3 = 2^2- Then the projection set V contains of 
Sil(Pi) = X2 

H2 = F3 = sn(A) = i 

G = 4X1 + + 6X1x1 - 3X| + 4X1 X| - 4X1 X| + 4Xf X2 + X| 

Note, that Zer:{G, M^) consists of two isolated points and and open curve. 
Our computations start and end with precision of 15 digits. The real in- 
tersection is computed in 0.305 seconds and consists of the isolated point 
(0,0,0) and an open curve (see Figure Ej). For the curve we computed the 
following three points. 

p = (1.91241422362700, -1.06499480841233,0.184566441477331) 
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Figure 7. A curve and an isolated point 

which corresponds to the Hft of the (non-isolated) critical point of Zer(G, M^), 
and 

(2.91241422362700, -2.54899069044757, 1.23313235073054), 
(2.91241422362700, -1.32006472767900, -0.443408797879122) 

which are sample points for the two branches ending and starting of p. 



Table 7. Experimental results for Example 15.181 



Size of Input 


Size of V 


Changes 


Precision 


Time 


1 


1 





15 


0.305 


4 


16 





26 


0.272 


8 


30 





40 


0.421 


12 


38 





48 


0.437 


17 


60 


7 


80 


6.119 


20 


72 





82 


1.215 


25 


92 


7 


112 


4.600 


28 


104 


1 


115 


1.900 


32 


118 





128 


1.716 


38 


134 


7 


154 


6.131 


41 


151 


6 


161 


5.303 


45 


165 


14 


194 


10.003 



2. COMPUTING THE REAL INTERSECTION OF QUADRATIC SURFACES 85 




Figure 8. Two intersecting lines with Sil(Pi) / 1 
Example 5.19 (a curve, Sil(Pi) / 1). Let be 

Pi = xl + xf-xl 

P2 = + XiX^ + — X3 + xf — X2 
P3 = X^ + X1X3 + X2X2 + X3 + xf — X2 

Then the projection set V contains of 

Sil(Pi) = xl-xl 

H2 = -1 + XI + X2 

= 1 — Xi + X2 

sn(K^) = xf-xl 

G = 1 

Our computations start and end with precision of 15 digits. The real inter- 
section is computed in 0.152 seconds and consists of two intersecting lines. 
We computed the following five points. The point 

p= (0,0,0) 

which corresponds to the lift of the critical point of Sil(Pi), and 

(-1,-1,0), (-1,1,0) and 
(1,-1,0), (1,-1,0) 

which are sample points for the two branches attached to the left and to the 
right of p. 
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Table 8. Experimental results for Example 15.191 



Size of Input 


Size of V 


Changes 


Precision 


Time 


1 


1 





15 


0.152 


4 


8 





18 


0.139 


8 


15 





25 


0.117 


11 


19 





29 


0.183 


15 


29 





39 


0.137 


19 


37 





47 


0.244 


23 


45 





55 


0.187 


27 


53 





63 


0.273 


31 


61 





71 


0.212 


35 


66 





76 


0.274 


39 


75 





85 


0.283 


43 


82 





92 


0.233 



Figure 9. One connected component 



Example 5.20 (one connected component). Let be 

Pi = ^2 - X3 + X1X3 + 5X2X3 + 2X| 

P2 = 6X| — 5X2X3 — X| + X1X2 — X1X3 + X3 

P3 = 6X| — 5X2X3 — X| + X1X2 — X1X3 + X3 
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Note, that P2 = Pj,- Then the projection set V contains of 

Sil(Pi) = ISX2 - I + 2X1 - Xl - 0^X2X1 - 2bXl 

H2 = H3 = Sn{Pi) = l 

G = -3X| - 8X^X1 -11X1 + 20X^Xf + 133X1 Xi 

+294X1 - X2XI + XfX2 + X2- X2X1 

Our computations start and end with precision of 15 digits. The real inter- 
section is computed in 0.529 seconds and consists of one connected compo- 
nent (see Figure [9]). 

2.7. Remark on Cubic Surfaces. We would like to remark that the 
algorithm presented in Chapter [2] has been extended to three cubic surfaces 
defined by the polynomials Ci, C2 and C3 in ]R[Xi, X2, X3]. Note that in this 
case the silhouette curve Zer(Sil(Ci), M^) contains all points (x, y) such that 
the polynomial Ci(x, y, X3) has a root z of multiplicity 2 or 3. Theorem 12.41 
implies that in the first case the polynomial sResi(Ci, 9Ci/(9X3)(x, y) 7^ 
whereas in the latter one sResi(Ci,(9Ci/5X3)(x,y) = 0. Moreover, one can 
also use a solution formula for cubic polynomials in one variable in order to 
lift a single point. 

Like in the case for quadrics, we can easily determine the missing adja- 
cency information while lifting the curve Zer(G, M^) using a simple combi- 
natorial type approach. 

Finally, this new algorithm has similarly implemented in the Computer 
Algebra System Maple (version 9.5) as well. The experimental results archived 
show a very good performance. We refer to |53j for more details. 
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